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ACCURATE evaluation OF TRANSIENTS ON 
DOUBLE CIRCUIT EH? TRANSMISSION LINES 

An overhead transmission line is subjected to overvoltage 
transients due to lightning, switching operations or faults. 

Among these, the overvoltages caused by the lightning is indepen- 
dent of the operating voltage of the line, whereas the over- 
voltages due to the switching operations or faults depend on the 
operating voltage of the line. Therefore, oflate as the trans- 
mission voltages are increasing, the studies of switching and 
fault transients have gained more importance. Accurate calcu- 
lation of these overvoltages permits reduction of safety margins 
and eventually reduces the cost of insulation of transmission 
lines and associated equipments. 

The various methods of calculating these overvoltage 
transients, described in the literature may be broadly classi- 
fied into the lumped parameter methods (including the Transient 
Network Analyser studies), time domain methods such as Lattice 
Diagram method, Uram and Miller's method and finite difference 
method and transform techniques. In this thesis, the modified 
Fourier transform method, first applied to the power system 
transient stMies by Mullineaux et al and later by Battisson 
et al and Wedepohl and Mohamad, has been made use of. Briefly 
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the Fourier transform method of calculating the transient 
involves solving the transformed differential equations of 
the transmission line for transformed voltage and current at 
any point on the line, satisfying the boundary conditions at 
the ends of the line hy the method of modal analysis and the 
transformation of the response to the time domain hy perfor- 
ming the truncated numerical iategration. The resistance 
and inductance matrices of a transmission line are in reality 
frequency dependent because of the variation of the conductor 
resistance with frequency due io the dkin effect phenomenon 
and the finite conductivity of the earth retiu’n path. Hence 
an exact stvidy of transients on transmission lines calls for 
the inclusion of the effect of frequency dependence of line 
parameters. The Fourier transform method, being a frequency 
domain method, has the special advantage that the frequency 
dependence of line parameters can be directly incorporated 
without involving almost any additional computational effort. 

The double circuit transmission line is a common 
feature of any power system. In the published literature, 
except for some TNA studies, the double circuit line transients 
have not been studied. Most of the double circuit lines have 
their conductors placed at mirror symmetrical positions with 
respect to a vertical plane. For a transposed double circuit 
line whose conductors are so symmetrically located, a 
frequency independent modal transformation is shown to exist 
in this thesis* For the untransposed double circuit line of 
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the abore structure, the eigen -value , eigen-vector analysis of 
two 3x5 matrices are to be performed at each frequency in the 
integration range to arrive at the modal transformation matrix. 

The objective of this thesis is to study the fault and 
energising transients on a typical double circuit line, using 
the modified Fourier transform method. 

The energisation and re-energisation overvoltages on a 
transmission line can be contained within the limit as low as 
1.5 p.u. by various control methods such as the m\ilti-step 
resistance closing, controlling the timing of the closure of 
the individual poles of the circuit breakers etc. When such 
effective methods of controlling the energisation and re- 
energisation transient overvoltages are employed, the over- 
voltages that arise due to fault initiation may become the 
limiting factor. Ninety percent or more of the faults 
experienced by a power system are single line to ground 
faults. In the event of a single line to ground fault, 
tremsient overvoltages of the order of 2 p.u. may be experien- 
ced by the sound phases. This calls for an accurate calcula- 
tion of transients due to single line to ground fault 
initiation. Various authors have studied this problem, 

Kimbark and legate have used the TNA to evaluate the fault 
initiating transient on a single circuit three phase line 
and have presented a Lattice Diagram approach for a theore- 
tical study of the problem with certain approximations. Later 
exhaustive analogue computer aided studies of the problem 
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have been reported. Bonnyubol et al have made use of the 
Laplace Transform method and evaluated the fault transient, 
employing the 'Residue Theorem* for the inverse Laplace 
transformation. Except for a few TRA studies, in all the 
above studies, only single circuit lineshave been consi- 
dered. Moreover the'- Lines have been assumed to be transposed 
and the effect of frequency dependence of line parameters has 
not been investigated. In this thesis, the modified Fourier 
transform method has been applied to study the transient due 
to the initiation of single line to grotmd fault at any inter- 
mediate point on a double circuit transmission line. The need 
for the modelling of a double; ;circuit line for fault transient 
studies has been established by comparing the fault transients 
on single and double circuit lines. A typical 400 KV double 
circuit line has been considered, talcing into account the 
mutual coupling between the two circuits. For the line, fed 
on both the ends by inductive sources, the mid span fault 
transients have been evaluated, assxining the line t o be both 
transposed and untransposed. The effect of frequency 
dependence of line parameters on the fault transient over- 
voltages for both the transposed and untransposed configura- 
tions of the line has been investigated. The fault location 
is changed to various off-centre points and the fault 
transients have been studied. The mid span fault has been 
fotmd to give rise to the maximum peak overvoltage. 

The energisation of a double circuit line involves 
closing of both the circuits. In practice, there is bound 
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to "be a time delay 'between the closing of the first and second 
circuits of a double circuit line. This delay is atleast a few 
seconds if not a few minutes. During this period, the transient 
due to the energisation of the first circuit would have died 
down and the first circuit would be in steady state. Thus 
essentially two transients are involved in the study of the 
energisation of a double circuit line (viz.) (i) the:- transient 
due to the closing of the first circuit, keeping the second 
circuit unenergised; (ii) the transient due to the closing of 
the second circuit when the first circuit is already in steady 
state. The thesis presents the calculation procedure and the 
results of the above mentioned studies on a double circuit line. 
A typical 400 K? double circuit line, open at the receiving end 
has been considered and the voltage waveforms at the open 
receiving end, subsequent to energisation from an inductive 
source have been calculated for both transposed and untrans- 
posed configurations of the line. The effect of the frequency 
dependence of line parameters has been incorporated in the 
studies by using the Carson's formulae for the calculation of 
the line parameter matrices. 



CHAPTER 1 


INTRODUCTION 

1.1 TRANSIENTS ON TRANSMISSION LINES 

The study of electrical transient phenomenon in 
power systems has gained great importance, especially with 
the increasing voltages of the transmission systems, The 
transient overvoltages are developed in a power system 
either due to internal or external causes. Lightning over- 
voltage is an example of the externally developed over- 
voltage and therefore the overvoltage caused by lightning 
is independent of the operating voltage of the ‘transmission 
system. Internally developed overvoltages are usually 
caused by the switching operations (either the opening or 
closing of the circuit breakers) and fault initiation or 
clearing. The overvoltages due to these depend on the 
operating voltage of the system and other system parameters. 
In earlier days when the system voltages were lower, the 
oveirvoltsges due to lightning were higher than those due to 
switching or faults. Nowadays, in high voltage systems, 
the overvoltages caused by switching or faults are higher 
than those due to lightning. Therefore, the basic insula- 
tion level (BIL) of a power system is being determined 
more by the behaviour of the system under switching and 
fault conditions than by its response to lightning surges. 
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Tile cost of insulation of equipment and transmission lines 
in a high voltage system is very high. When the BIl cannot 
be determined accurately, the factor of safety has to be 
rather high. However, accurate determination of the BIL 
permits the reduction of this factor of safety and thereby 
the cost of insulation can also be made an optimum. Further, 
once the magnitude of the overvoltages that can arise in 
a system is known, methods of controlling them within limits 
can be developed. An accurate analysis requires an exact 
modelling of the system components and the physical pheno- 
mena involved in the study and similarly an accurate method 
of solution has to be employed. 

When a switching operation or fault takes place, 
the elements of the power system are subjected to voltages 
and CTirrents, having a wide frequency range, which may 
extend from 50 Hz to 100 KHz. The resistaiace and induc- 
tance matrices of a transmission line are, in reality, 
frequency dependent due to the skin effect phenomenon and 
the finite cohductivity of the earth return path. Many 
of the EHV transmission lines are untransposed and even 
if a line is transposed, an ideal transposition is not 
practically realizable. Hence an exact study of transients 
on transmission line calls for the inclusion of the effect 
of frequency dependence of line parameters and handling of 
the untransposed lines. 
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There are various methods of calculating these 
overvoltage transients on a transmission network, avail- 
able in the literature. These may be broadly classified 
into lumped parameter method, lattice diagiram method, 

Uram and Miller's method, Fourier transform method, 

Dommel's method and finite difference method. Brief 
description of these methods and their salient features 
are presented in the following sections. 

1.2 LUMEBD PARAMETER METHOD 

In the lumped parameter method, the transmission 
line is represented by a number of lumped T or tc sections 
and the differential equations of the resulting equivalent 
network are solved either with a transient network analyser 
or a digital computer. 

The favourable features of the method are its 
simplicity and its ability to handle imtransposed lines 
and nonlinear elements directly. 

On the other hand, for achieving good accuracy in 
the results, a large nmber of line sections are required 
in this method and this correspondingly increases the 
storage requirements. The frequency dependence of line 
parameters cannot be incorporated in the method. Further, 
the computations are to be carried out in stages right 
from the initial state. In other words, the response at 
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any particular instant cannot be directly found out without 
going through the calculations in stages from the initial 
state. 

1.3 LATTICE DIA&RAM METHOD 

The Lattice diagram was first introduced by 
Bewley [l]. The application of this method to large 
scale single and three phase systems has been done by 
Barthold and Garter [2], Bickford and Doepal [3] and 
others. In this method, lines and cables are specified 
by their surge impedances and surge travel times and the 
reflected and refracted voltages and currents at junctions 
and terminations are calculated by the use of reflection 
and refraction coefficients. To reduce the computational 
effort, the l\imped parameters such as the transformers, 
shunt reactors etc. are represented as transmission line 
stubs of appropriate reactance values with the terminals 
kept open or shorted. 

The advantages of this method are that it is 
simple, takes less computational effort and can handle 
large systems. 

On the contrary, in this method, the frequency 
dependence of line parameters cannot be incorporated. 

Line losses can be taken into account only approximately 
by an attenuation factor. Representation of lumped para- 
meter elements by stub lines leads to inaccuracy if proper 
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values of surge impedance and travel time for these stubs 
are not chosen. A lot of data such as the reflection 
coefficient array, travel time array, the magnitudes of 
the forward and backward waves are to be stored. Hence 
the memory requirement is quite high. 

1.4 UIIAM AND MILLER'S METHOD 

The Uram and Miller's method [4,5] is a time domain 
method based on Laplace transform technique. In this 
method, by making certain assumptions, an incremental 
solution of the transmission line equations is made 
possible. Here the voltage and current at each instant 
of time and at any point on the transmission line, subse- 
quent to a switching operation, can be expressed in terms 
of the values of some four functions at that instant 
(which can be interpreted as the forward and backward 
waves). Among these four functions, two are determined 
by the terminal conditions existing at the ends of the 
line and the other two are just the attenuated delayed 
functions of these two fuactions. For the case of a line, 
which was unenergised before the switching operation, the 
delayed functions take zero values initially till the time 
equal to the travel time of the line and this facilitates 
the starting of an incremental solution of the problem. 

For the three phase case, the functions become vectors of 
order 3, the components of which correspond to the modal 
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quantities. The method basically uses Laplace transform 
technique. The transformed transmission line equation is 
decoupled into modal components with the help of a constant 
modal transformation matrix. Such a transf oimation exists 
only for the transposed line and hence the applicability of 
the method is restricted to the transposed lines only. With 
the other assumption that the resistance offered by the 
line for the various modes is much less, when compared to 
the inductive reactance offered by the line for the various 
modes at all frequencies, the inversion of the response in 
the Laplace domain to time domain by analytical means is 
made possible. 

The favourable features of the method are that it 
involves less computational effort and takes less computer 
time. 

But this method cannot take into account the effect 
of frequency dependence of line parameters directly. The 
untransposed lines cannot be handled by the method. Compu- 
tations are to be carried out in stages, starting from the 
initial state, Wonlinearities cannot be handled directly, 

1.5 COURIER TRANSFORM METHOD 

The Fourier transformi method has been applied for 
studying the transients on transmission lines by Day, 
Mullinaux et al [6-8], Battisson et al [9] and Wedepohl 
and Mohamed [lO], The method involves solving the Fourier 
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transformed differential equations of the transmission 
line for transformed voltage and current at any point on 
the line, satisfying the boundary conditions at the ends 
of the line by modal analysis and the transformation of 
the response to time domain by numerically performing the 
inverse integral, A modified version of the Fourier 
transform, which includes a negative exponential term in 
the basic transform pair relationships makes the integral 
over the infinite range numerically stable. The infinite 
integral is truncated at a finite large value for the 
purpose of performing the numerical integration and the 
resulting xinwanted oscillation (G-ibb’s oscillation) in 
the process is avoided by introducing a factor to the 
integrand. As the method involves calculation of the 
frequency response in the whole frequency range, the 
additional calculation of the modal transformation matrix 
at each frequency to handle the untransposed line case 
and the calculation of the line parameter values at each 
frequency to be used in the transients computations are 
s t r a ight f o rwar d . 

This method does not involve any assumption and 
represents all the components in a detailed manner. 
Incorporation of the frequency dependence of line para- 
meters in the method is direct and does not involve any 
additional computations but for the calculation of just 
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the parameters at each frequency. Calculation of the 
solution just at any desired instant or period can be 
done and does not have to be carried out in stages right 
from the initial conditions. Handling of lumped parameter 
elements is very simple. There is no difficulty in 
studying the untransposed line cases by doing the eigen- 
value, eigenvector analysis at each frequency to arrive 
at the modal transformation matrix at that frequency. 

1.6 DOMMEI’S METHOD 

The Dommel's solution procedure [11,12] adopts the 
method of characteristics [15>14] for the transmission 
line and the trapezoidal rule of integration for the 
lumped parameter elements. According to the method of 
characteristics, for an imaginary observer, moving along 
the transmission line in the forward direction with a 
velocity equal to that of the travelling waves initiated 
on the line, following a disturbance, the expression 
(v + z i) remains constant, where v and i are the voltage 
and current respectively at any point on the line and z is 
the surge impedance of the line. Also for an observer, 
moving along the line with the same velocity in the back- 
ward direction, the expression (v - z i) remains constant. 
Forcing these conditions, the current and voltage at both 
the ends of the line can be expressed in terms of the 
surge impedance of the line and the past values of the 
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voltages and currents at these terminals, prior to a 
period equal to the travel time of the line. The appli- 
cation of the trapezoidal rule of integration for the 
lumped parameters facilitates expressing the currents and 
voltages at the terminals of these elements at each instant 
of time in terms of the values of these currents and 
voltages an j.nterval earlier and in terms of the values 
of the parameters. Once each component of the network is 
replaced by its equivalent circuit described above, the 
nodal equations for the whole network at each instant of 
time can be written and they can be solved for the unknown 
node voltages, knowing the past history of the system over 
a time span equal to the travel time for the transmission 
line and just the previous step value for lumped parameter 
elements. The line losses can be represented by lumped 
resistances at the ends and at the middle of the line, 

For solving the nodal equations of large scale networks, 
the application of the well developed optimally ordered 
tri angular! zat ion and the sparsity techniques reduces the 
storage requirement considerably. 

This method is well suited for large scale system 
problons especially because of the availability of the 
well developed sparsity techniques. It can handle non- 
linearities. The method is quite general and hence the 
development of a general purpose program is possible. 
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On the other hand, there is no direct way of incor- 
porating the frequency dependence of line parameters in this 
method. Computations are to he carried out in stages, 
starting from the initial state. As the method represents 
the lossy nature of the line by lumped parameter elements, 
to that extent the modelling is not accurate. 

1.7 FINITE DIFFERENCE METHOD 

The finite difference method has been applied for 
solving the electrical transient problems by Stafford 
et al [15] and Raghavan and Sastry [16,17]. In this method 
for determining the voltage and current on the transmission 
line, which is a function of both distance and time, the 
partial differential equations of the line are converted to 
difference equations, adopting a central, forward or back- 
ward difference scheme. The distance— time space of the - 
transmission line is divided into a number of lattice 
points and the difference equations are solved simultane- 
ously step by step to arrive at the voltage and current at 
each lattice point, satisfying the boundary conditions at 
the ends of the line. 

The nonlinearities and untransposed lines can be 
handled by this method. 

However, proper choice of time step and distance 
step is very crucial for ensuring nximerical stability. The 
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method requires large storage, Frequency dependence of 
line parameters cannot be directly incorporated. 

1.8 OBJECTIVES AED OUTLINE OP THE THESIS 

The objective of this thesis is mainly to develop 
and extend the modified Fourier transform method for 
solving the electrical transients on double circuit trans- 
mission lines, A double circuit line is a part and parcel 
of any power system. In the existing literature, the 
electrical transients on double circuit lines have been 
studied mostly using the Transient Network Analyser. This 
thesis presents the development of a detailed model for a 
double circuit line and analyses the transients on them 
due to fault initiation and line energisation by the modi- 
fied Fourier transform method. The modified Fourier 
transform method has been rated as one of the more accurate 
methods for solving these problems as it can be applied to 
a detailed transmission line model with ease. The distri- 
buted nature of all the parameters of the line is preserved 
in the analysis. The effect of frequency dependence of 
line parameters can be accounted for easily in the method 
almost with no additional computational effort. The case 
of the untransposed line can also be handled by this 
method. 

Also the relevant theory for extending the Uram 
and Miller’s method for studying the electrical transients 



12 


on transposed double circuit transmission lines is deve- 
loped and presented in the thesis. As this method assumes 

1 

a simpler model for the transmission line, it takes 
relatively less computer time and gives fairly accurate 
results. However, it can be used only for a transposed 
line and without incorporating the effects of frequency 
dependence of line parameters. 

In this chapter a general introduction to the 
study of electrical transients on transmission systems 
has been given together with a brief description, and the 
salient features of the various methods that have been 
reported in the literature for evaluating these transients. 

Chapter 2 describes the theory of the modified 
Pourier transform method as applied to the study of elec- 
trical transients on three phase power systems. The concepts 
of matrix functions and modal decomposition, proposed by 
Wedepohl and Mohamed [10], are made use of in the solution 
procedure. A typical 400 KV single circuit line, open at 
the receiving end is considered and the transient voltage 
waveform at the receiving end arising due to its energisa- 
tion is computed, A comparison of voltage waveforms at 
the open receiving end for the transposed and untransposed 
configurations of the line is done. 

A mathematical model of a double circiiit transmission 
line, suitable for electrical transient studies is developed 
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in Chapter 5. It is shown in this chapter that for a 
transposed double circuit line of mirror symmetric 
conductor configuration with respect to a vertical plane, 
there exists a constant modal transformation matrix, which 
holds good at each frequency of the integration range. 

An EHV transmission line is subjected to high 
overvoltages due to its energisation and re-energisation 
on trapped charge. Effective ways of restricting these 
overvoltages to low values, such as multistep resistance 
switching, controlling the timing of the closure of the 
individual poles of the circuit breaker etc. have been 
developed and are implemented in practice to limit these 
overvoltages to less than 2 p.u. When such methods are 
employed for controlling switching overvoltages, the over- 
voltages that arise due to fault initiation and fault 
clearing may dominate the scene over the switching over- 
voltages. Because of this, the study of fault transients 
on transmission lines is of importance. In Chapter 4, the 
application of the modified Fourier transform method for 
evaluating the overvoltage transients, arising on the sound 
phases of a transmission line due to a single phase to 
ground fault at any intermediate point on the line is 
described. The mathematical model, developed in Chapter 3 
for a double circuit line is used to study the transient 
overvoltages on the sound phases due to a single phase to 



14 


ground fault on one circuit of a typical 400 KV double 
circuit line. A comparison of fault overvoltages for the 
transposed and \mtransposed configurations of the line is 
made. The effect of frequency dependence of line parameters 
on the fault initiating tiansient is studied for both the 
transposed and untransposed double circuit lines. The 
fault transient is evaluated for various intermediate 
fault locations. 

The modified Fourier transform method is applied 
to compute the energisation transients on a double circuit 
line in Chapter The model, developed in Chapter 3 for 
the double circuit line is used in these studies. The 
energisation of the first circuit with the second circuit 
remaining unenergised as well as the energisation of the 
second circuit when the first circuit is operating in 
steady state are considered. In these studies, it is 
assumed that all the three poles of the circuit breaker 
of each circuit get closed simiiLtaneously . This simpli- 
fies the analysis and helps in getting an idea of the 
magnitudes of the overvoltages involved, A 400 KV line 
example is considered and the energisation transient at 
its open receiving end is computed for both the transposed 
and untransposed configurations of the line. 

The extension of the Uram and Miller’s method of 
solving the electrical timisients for a double circuit 
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line is done in Chapter 6, The existence of a constant 
modal transformation matrix, independent of frequency, 
for a transposed double circuit line, as shown in 
Chapter 3» has made the application of this method for 
this problem possible. This being a time domain method, 
it takes much less computer time but can handle only 
transposed lines and the frequency dependence of line 
parameters cannot be directly incorporated. The method 
is applied to the same 400 KV double circuit line example 
as of Chapter 5. The results, obtained for the example 
for both the first and second circuit energisation are 
presented. 

In Chapter 7, the effect of the nonsimultaneous 
closing of the circuit breaker poles on the energisation 
transients is investigated. The theory of the modified 
Fourier transform method is extended for analysing the 
energisation transients on double circuit lines, 
incorporating the sequential closing of brea,ker poles. 

The same 400 KSf double circuit line example as of 
Chapters 5 and 6 is considered and the energising tran- 
sients are evaluated for both the transposed and untrans- 
posed configurations of the line. The effect of frequency 
dependence of line parameters on the transients is 
investigated. 
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In. the concluding chapter 8, the results of the 
various chapters of the thesis are revievred, conclusions 
are drawn and problems for further research are outlined. 



CHAPTER 2 


MODIFIED FOURIER TRANSFORM METHOD 

2.1 INTRODUCTION 

The theory of the modified Fourier transform method 
and its application for solving the electrical transient 
problems [9flO] have been described in this chapter. This 
method has been adopted in this work with a view to repre- 
sent the transmission line and other components in a 
detailed manner so that the transients can be evaluated 
more accurately. It is easy to incorporate the effect 
of frequency dependence of line parameters in this method 
of transients evaluation. The method takes all the para- 
meters of the transmission line to be distributed and can 
handle untransposed lines. To illustrate the method, an 
example of a single circuit line has been considered and 
the results of the energisation study are presented, 

2.2 THEORY OF THE METHOD 

The solution of the transients on multiconductor 
transmission lines is complicated by the fact that the 
partial differential equations describing the behaviour 
of the system are 4-dimensional (3-dimenaions in space 
and the time dimension). Ry making the assumption that 
only plane-wave propagation takes place, it is possible 
to separate the two space variables mutually at right 
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angles to the direction of propagation and to express them 
in terms of equivalent impedances, as described by Oar son 
[18], This leaves space variable in the direction of 
propagation and the time variable alone in the analysis. 
With this assumption, the voltage and current at any point 
on a transmission line are determined by the partial 
differential equations 


9 v(x,t ) 

8 X 


= R. i(x,t) + L 




3t 


( 2 . 1 ) 


_ ( 2 . 2 ) 

■where v(x,t), i(x,t) are the voltage and current respec- 
tively at the point, x units away from the end at time 't' 

R is the resistance/unit length of the line 

1 is the inductance/unit length of the line 

G- is the conductance to ground/unit length of the line 

and C is the capacitance to gromd/unit length of the line. 

For a transmission line, usually G is negligibly 
small and hence (2,2) takes the form 

= c (2.3) 

When equations (2,1) and (2.3) are Fourier trans- 
formed, they become full differential equations in trans- 
formed voltage and current. 

The Fourier transform pair is defined by the 


expressions 
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I'(m) 

= / f (t) dt 

0 

(2.4) 

f(t) 

= -^ / F(m) do) 

(2.5) 


Equations (2.1) and (2.3), when Fourier transformed take 
the form 

_ IZkiLiSli = [H + ouL] I{x,m) (2.6) 

- V(k,Ii)) (2.7) 

Differentiating equation (2,6) with respect to x and 
substituting for from equation (2.7), we get 

^ = [r + jojL] 3 wC 7(x,(o) (2.8) 

dx^ 

Solution of the above equation can be written as 

V(x,u)) = Cosh yx.‘A + Sinh yx*® (2.9) 

where y = ([R + jwlljuC)^ 

and A and B are arbitrary constants, which are determined 
by the terminal conditions, existing at the ends of the 
line. 

From equation (2,6), the expression for trans- 
formed current can be obtained as 

I(x, w) = -[ ||^- |;j3^[Co8h yx-B + Sinh yx-A] (2,10) 
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Once the boundary conditions are known in the 
frequency domain, the arbitrary constants A and B can be 
evaluated and the transformed voltage and current at any 
point on the line can be found out with the help of 
equations (2. 9) and (2.10). As this response in frequency 
domain is a complicated expression involving hyperbolic 
functions, it can be transformed back to time domain only 
by numerically evaluating the inverse integral of the 
form (2.5). The above mentioned conventional Fourier 
transform technique is not suitable for the evaluation of 
this integral numerically, as the integrand in this foim 
is not quite stable numerically. A negative exponential 
term, when introduced to the integrand makes it numerically 
stable and hence a modified version of the transform is 
preferred for the calculations. 

The modified Fourier transform pair is given by 
the expressions 

V(x,a)) = /v(x,t) dt (2.11) 

0 

v(x,t) = ^ / V(x,w ) dm (2,12) 

where a is a positive constant, sufficiently large to 
ensure that the above integrals converge. 

When the modified Fourier transform is applied 
instead of the conventional Fourier transform. 
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jquations (2.6) and (2.7) get modified to the form 


ctV(x«a>l = [j^ + (a + 0 (o)L] l(x,w) 
dx 


^ <tl(x, ml ^ (a + jw)C V(x,w) 
dx 

Putting these equations together, we get 


d^V(x, ml _ 2, y V(x,m) 


( 2 . 15 ) 


( 2 . 14 ) 


( 2 . 15 ) 


where 


Z = R + (a + jw)! 
Y = (a + ow)C. 


The general solution for transformed voltage and 
current at any point on the line can now he written as 

V(x,w) = Gosh **■ Sinh 

I(x,u) = -(T/Z)^ [Cosh yx-B + Sinh yx-*.] (2.17) 

where now y ={ [E + (a + im)I.]Ca + • 

The boundary conditions, when substituted in 
equations (2.16) and (2.17) determine the constants A 
and B. The final step is to evaluate the integral (2.12) 
to arrive at the response in time domain. 

Por evaluating the above integral numerically, the 
infinite range is, in practice, truncated at a large finite 
value, say S . To effect economy in computation, a half 
range foxt, of the Fourier Integral can be used and in that 
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case the integral to he evaluated takes the form 
p at 

v(x,t) = — I Real Part[V(x, w) JCosiot dm 
^ o 

The truncation of the infinite range of the integral 

at a finite value has been found to give rise to unwanted 

oscillations» known as G-ibb*s oscillation in the solution. 

To eliminate these oscillations, a wei^ting function known 

as sigma factor is introduced to the integrand. With the 

introduction of this factor, the integral reads as 

p .at Q Sin( ^ 

v(x,t) = — / Real Part [V(x, ta) ] Cos wt — (2.18) 

° ) 

A set of suitable parameters (J2,Aw and a) to achieve 
accuracy and economy in computation, while evaluating 
numerically the above integral has been reported [6-8]. 

Accuracy in representation of the rise time of the 
response is determined, in main, by the choice of . For 
a good representation, it is significant to choose so that 
the calculated rise time is less than 10 percent of the 
travel time of the transmission line. 

i . e . T < t/ 10 

It has been found, in practice, that the inaccuracy in the 
rise time is principally in the range Si T > 4. Bo the 
above criterion implies choosing Si so that 

n T > 4 

i.e. Si > 40/t 
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It has been confirmed by experience that the arjgular 
frequency interval Am = n/500 is a good choice. 

While the paarameter 'a* can be used to 'smoothen* 
the term V(a,w) in the integrand, the ' Ooswt' term oscil- 
lates with a period in m of Zn/t, For tahing into account 
this oscillation, introduced to the integrand accurately, 
it is desirable to choose Aoj, so that at least 4 points 
occur in each quarter cycle 
i .e . 16 Am < 27i/t 

In other words, the period over which the response 
is calculated should be restricted for a particular value 
of A(o chosen so that 

period of calculation < 0.4/ Aw 


2.3 THREE IHASE SYSTEMS - MODAl ARALTSIS 

The theory involved in applying the modified Fourier 
transform method for solving the transients on single phase 
systems has been presented in the previous section. When 
this is extended to a 3 phase system, the voltages and 
currents become vectors and the line parameters become 
matrices. 


Now the matrix differential equation, which has to 
be solved to arrive at the transformed voltage vector at 
any point on the line takes the form 
d2V(x,a)) 


dx' 


[P] Kx, m) 


( 2 . 1 ^) 
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where [P] = [z][Y] 

[Z] = [R] + (a + ja3)[l] 

and [Y] = (a + j (a)[C] 

For solving the coupled matrix differential equation 
(2.19), the concept of matrix function is to he made use of- 
By the concept of matrix function, if a mathematical 
operation is to be performed on a square matrix, the resul- 
ting matrix after the operation can be obtained by doing 
the operation individually on the eigenvalues (modal values) 
of the matrix and by pr e-mult i plying and post -multi plying 
the diagonal matrix, having these operated eigenvalues as 
diagonal elements, by the eigenvector (or modal transforma- 
tion) matrix and its inverse respectively. 

Symbo li cally , 

f([A]) = [Q] f([X]) [Q]-l 

Using this property, the solution of equation (2.19) 
can be written as 

V(x,a)) = [Cosh'I'x] A+ [Sihh^ x] B (2.20) 

where [t|)] = 

[Cosh ^x] = [Q][Gosh 

[Sinh i|»x] = [Q][Sinh yx] [q]"^ 

[y] is the diagonal matrix, whose diagonal elements are 
the square roots of the eigenvalues of matrix [P] = [z][Y] 
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[Q] is the modal transformation matrix, which makes 
[Q]’’^CI’][Q] diagonal 

and A and B are arbitrary constant -v-ectors, which can be 
evaluated by the terminal conditions of the line. 

Ihe general current vector solution can now be 
determined as 

I(x,u)) = ]{[Cosh t|ix] B + [Sinh ij^x] A } (2.21) 

When the boundary conditions are forced into 
equations (2.20) and (2,21), the constant vectors A and B 
are determined and hence the transformed voltage and 
current vectors at any point on the line is completely 
determined. Then the final step of evaluating the inverse 
integral of each of the phase components of the transformed 
response separately for the 3 phase case yields the phase 
voltages and currents in time domain. 

In the above analysis, if the transmission line is 
transposed, the computation of the eigenvalues y’s 
eigenvector [Q] becomes unnecessary, as for this case a 
constant modal transformation, which holds good at each 
frequency in the integration range exists, for the trans- 
posed line, the line parameter matrices and hence the 
matrix [P] have a special structure, in which all the 
diagonal elements are equal and all the off-diagonal 
elements are equal. GJhough the values of these diagonal 
and off-diagonal elements vary with frequency, still at 
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each frequency this special structure of equality of 
diagonal and off-diagonal elements is preserved. The eigen- 
values of such a matrix [P] are the one corresponding to the 
ground or zero sequence mode, given hy (diagonal element + 

2 X off-diagonal element) and the other two (repeated ones) 
corresponding to the aerial modes given hy (diagonal 
element - off-diagonal element). An eigenvector matrix for 
the matrix [P] of such a structure [4] is given hy 



( 2 . 22 ) 


When the line is untransposed, the matrix [Pjdoes 
not possess this special structure. Hence for this case 
the eigenvalue, eigenvector analysis has to he performed 
at each frequency to arrive at the matrices [Q] and [y]» 
which hold good only at that frequency. 


2 ,4 mathematical model for bnergisatioh study 

In this section, the development of an algorithm 
for computing the transient following the energisation of 
a single circuit, three phase line [ 10 ], applying the 
theory of the previous two sections is described. The 
algorithm assumes simultaneous closing of the circuit 
breaker poles during energisation. 
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Supposing the transformed voltages at the sending 
and receiving ends of a transmission line are known, forcing 
of these conditions in the general voltage vector solution 
expression (2,20) yields 

A = Is (2.23) 

[Gosh4»Jl]Vg + [Sinhi};Jl]B = 
where 5. is the length of the line. 


i.e. B = [Cosech ipA - [Goth (2.24) 

Prom (2.21), (2.23) and (2.24), the transforms of currents 
at the terminals of the line can be expressed in terms of 
the terminal transformed voltages as follows: 


I 

~s 


[Ylll 

[X^] 


Y ~ 

— s 



j 1 1 

1 

1 1 

[Xi^ 




where t^ll^ = t(>][Coth 

and [Tj^] = -[Z]"^[ 'l’][Coseoh.i|»l] 


(2.25) 


In the above equation, the matrix relating the 
terminal voltages to the terminal currents can be inter- 
preted as the short circuit admittance parameter matrix of 
the line, viewed as a two port network. 


Supposing a shunt element is added at one of the 
terminals of the transmission line, then the short circuit 
admittance matrix of the resulting overall set up can be 
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obtained by just addiiig to the self admittance term, 
corresponding to the particalar terminal of the transmis- 
sion line, the admittance valae of the added shunt element. 


Por simulating the energisation of the line from a 
constant voltage source having a source impedance, conversion 
of the voltage source with the series impedance to its 
equivalent current source with the appropriate shunt admit- 
tance facilitates the application of the modified short 
circuit admittance parameter matrix (obtained as described 
above) in the analysis. 

Figure 2 .1.1 and 2.1.2 give the system set up 
with voltage source and with the equivalent current source 
respectively. 


Referring to Figure 2.1.2 when the shunt elements 
corresponding to the source admittances are added at the 
sending end terminals, the short circuit admittance para- 
meter of the system takes the form 


where 


(2.26) 


l/[Rg+( 

0 l/[Rg-t-(a+im)l.g] 0 

0 0 V[Rg+(a+jai)Lg] 


(2.27) 



— c> 


S 2.1,1s SYS? ‘IK DIAGRAM WITH SODRCB ERPHBSENTATIOIT 


[Rg+(a+3w)L«3 
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Now the desired response is calculated for the 
injected current vector at the sending end 


I 

“S 


E/[Rg+(a+3w)lg] 

S(,/[Rg+(a+ja))Lg] 


(2.28) 


where E , E and E_ are the modified Fourier transform of 
A j3 C 

the source voltages of the three phases, which are analy- 
tically known. 


Considering the case of a line, open at the receiving 
end, equation ( 2 . 25 ), when the boundary conditions are sub- 
stituted take the form 




-s 


g 



([Yu] + [Y 3 ]) 




i — 1 




Is 

[Yu] 


A 


(2.29) 


Prom the above equation, expressions for and can be 
obtained in terms of I as follows ; 


Vf, = is (2.51) 


2.5 NUMERICAI. EXAMPLE 

The above algorithm is applied to study the 
transient following the energisation of a typical 400KV 
single circuit, three phase line, open at the receiving 
end and the results of the study are presented in this 
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section. The source has been assumed to be purely 
inductive. The effect of the presence of earth wire is 
included by mating the necessary correction to the 
phase parameter matrices of the line. The line parameters 
and other system data are given in Appendix A. 

The phase voltage waveforms at the open receiving 
end, following the energisation for the transposed and 
untransposed configurations of the line are compared in 
Figure 2.2. All the three poles of the circuit breaker 
are assumed to close simultaneously at the instant when 
the source voltage of phase- A is at its peak value. A 
maximum peak of -2,45 p.u. is observed to occur on phase- A 
at 9.2 m.secs. for the case of the transposed line, whereas 
the maximum peak for the case of the untransposed line is 
2,22 p.u. on phase- C at 13.26 m.secs. 

2 .6 CONCLUSIONS 

The theory of modified Fourier transform method 
for studying the electrical transients on transmission 
lines has been presented in this chapter, A typical 400KV, 
single circuit, three phase line example has been consi- 
dered and the transients following its energisation from 
an inductive source has been evaluated for both the cases 
of transposed and untransposed configurations of the line. 
The extension of fhe theory for evaluating the energisation 
and fault initiating transients on double circuit transmis- 
sion lines is described in the following chapters. 


CHAPTER 5 


MODELLING OP A DOUBLE CIRCUIT LINE 
FOR transient STUDIES 


3.1 INTRODUCTION 

In this chapter, a model of a double circuit 
transmission line, suitable for transient studies is 
developed. This model is used in the subsequent chapters 
for analysing transients due to energisation and fault 
initiation. The double circuit line is a part and parcel 
of any power system. So far, in the literature, except for 
a few studies using Transient Network Analyser (TNA) [19] » 
the double circuit line transients have not been reported. 

The study of transients on a transmission line 
involves solving of the polyphase x^ave equations. The 
first step in the solution procedure is to decompose the 
equation to modal components. The transformation matrix 
for this purpose is to be obtained by performing the eigen- 
value, eigenvector analysis of the line propagation coeffi- 
cient matrix. 

In any practical power system, most of the double 
circuit lines have got a mirror symmetric conductor 
configuration with respect to a vertical plane. This 
gives most economical tower configuration. For such a 
special conductor configuration, a saving in the compu- 
tation of the eigenvalues and eigenvector can be 
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aficbhio iwed [20]. In "the <!a_se of a:.ii ifie^ly transposed 
fdoaubltls circuit line orf nirr^orr ^jr:miettrlc conductor confi~ 
^ima’jitiuon, a corstan-t ao da-l tiranis -foorma-t ion matrix, which. 
lilioIMiJs good at aLl fareqaeaciLes ca_D b® shown to exist. 

S'S2 1 DDAL TRAlTSPOHIJlOO I HOH .i DOlBBLiE CIRCUIT 
■ eiRSMISSION IIIB 

The line paraU'et er“ aoa"tr;i«eB af a double circuit 
flrs'aamsaaission Line, ’wt-ose ooncdu_ct=orr oca nfiguration is mirror 
3 Bji'iiaiieWtric with respect to a "re: rteicuall plane are synuaetrical- 
Ilhue isatrix [T] of the llne^ ■•whJ.ccli ias given by [z][Y] vilIL, 
t th-ifiirrerfore also be 3 yIaJne^tl^ica 1. ?L-Jir~tker the matrix [P] 
oor'’tdthMe line parameter mfcric e», ntfaern partitioned, take 


km Oecial form 


II 

1 1 

[M^] 



ijM^] 

[M,2 


IitfutHL-s case of htis a .37“ toe mri' .tfcen as 

3 C 
3 I 

(Itiiemj® [P ] and [E-l] are natiriccesa ofr order 5- These 
sputaiiwtrices [T ] and [X] ar*e s; 3 'iuiiestx'ical. However, 
ttlifiisirr diagonal elements ar“e iiB3e*<jiLiaH< Similarly , their 
otff->®iagonal elementa ar e a J-sao * uoDeena-l . 
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The eigenvalues of the matrix [P] are obtained by 
solving the poUynomial equation 

det([P] ~ [U2]) = 0 (3.2) 

for the six values i = 1,2,..., 6, where [U2] is unit 
matrix of order 6. 

Since the determinant of a matrix product is equal 
to the product of the determinants of the component matrices, 
taken in any order, if we choose a particular regular matrix 



[u] 

-[u] 

[K] = 


[U] 


where [TJ] is unit matrix of order 3» the equation for 
finding the eigenvalues can be written as 

det{[K:]”^([P] - Xi[U2]) [K]} = 0 (3-4) 

It can be verified that the matrix product 

[ 0 ] 


^0 0 
It 

where [O] bp 0 0 

^ 0 0 

X^^; i = 1,2,3 are the eigenvalues of the matrixCC^g^l+El^j] ) 




;(Cpj+[Pt]-x,i[u:j 


[K]“\[P]-Xl[D2])M = 
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and i = 1,2,3 are the eigenvalues of the matrix 

([Pa]-[Pt]). 

Therefore, equation (3*4) becomes 

([Pj+[P,]-X3i[U]) [0] ^ 

det - 0 

COJ ([Pj-[Pt]-^4^ 

i.e. det ([P^]+[P^]-X^.[U]) det ([P^]-[P^]-X^^[tr]) =0 

(3.5) 

Hence the eigenvalues of the matrix [P] can be 
obtained by solving the equations 

det ([P^]+[P,,]-Xgj^[U]) = 0 for j 1=1,2, 3 

and det ([Pg]-[Pj,]-X^^[U]) = 0 for X^^^ ; 1=1, 2, 3 

separately. 

The eigenvector matrix [Q] of the matrix [P] is 
to be obtained by solving the matrix equation 

[P][Q] - [X][Q] (3.6) 

where [X] is the diagonal eigenvalue matrix of matrix [P]. 

Equation (3*6) can be manipulated to the form 

[Ertp][E][K]-^[Q] = [K]-^[X][E][k]-^[Q] 

where the matrices [K] and are the same as described 

earlier. 
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¥e have seen that 


[K]-1{[P]-[X])[K] 


([Pj+[Pt]-[X3]) [0] 

[0] ([PaWf^l-M) 


■where [Xg] is the diagonal eigenvalue matrix of the 
matrix ( [ ^ diagonal eigenvalue 

matrix of the matrix ([P„]-[E^]). 

D 

Let [K]“^ [Q] be denoted by [Q’ ]. Then 


([Pj+CPfeKx^]) 

[0] ' 


[Qii] 

[ Q5l2 3 

[0] 

([Pj-CP^K^a])- 


[Q 21 ] 



, 

[0] 

[0] 

tBS 

£0] 



(3.7) 


where [Q^], and are submatrices, 

obtained after partitioning of matrix [Q'j for conformable 
muj-tiplication in the above equation. 


Prom equation (3.7), it is evident that 

if det([Pg^]+[P^]-[Xg]) = 0 

and det([P^]t*[P^]-[X^]) ^ 0 , 

then *= 

and is obtained by solving the matrix equation 

( [p^]+C5jj]-C Xg]) = [o]. 
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SimilaJ-ly if det( [P^]-[P^3-[ X^] ) = 0 

and det( [P^]+[P^]-[ Xg]) 4 0 , 
then [Q^] = [O] 

and obtained by solving the matrix equation 

([Pa]-[P^]-[x^])[Q^2^ = 

Thus [Q*] takes the form 



The above form of the matrix [Q' ] is very useful. 

It contains only the self term [0223* mutual 

terms [0^3 and [Q2 i 3 are zero. This fact is of particular 
importance in the present calculations. Therefore » 


[q3 = [k3[q'3 = 


[u3 

-[IT3 

[Qii3 

[03 ^ 

[u3 

[tJ3 

; [03 

[0^2^ 


3 - 

3 

The above matrix [q 3 again contains only the self 
term [QJj_ 3 and [ 0 ^ 2 ^ mutual terms [Q^3 and [Q^3^3 

are not present. Therefore, the problem of solving the 
transient on a double circuit line can be reduced to the 
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solution of two single circuit lines, decoupled from each 
other. Thus, all operations and manipulations can be carried 
out independently on each circuit as if the other does not 
exist. Only when the composite parameters are required, the 
independent values are combined by means of matrix [K]. 


3.2.1 Transposed line; 

For an ideally transposed line, the submatrix [P ] 

takes the particular form, in which all the diagonal elements 

are equal and all the off-diagonal elements are also equal. 

In addition, all the elements of the submatrix are also 

equal. Hence the matrices ([P ]+[P. ]) and ([P ]“[P, ]) will 

a b a b. 

also have the structure, in which all the diagonal elements 
are equal and all the off-diagonal elements are equal. The 
diagonal and off-diagonal elements of these matrices 
([P^]+[P^]) and ([Pg^]-[P-jj] ) (which are dependent upon the line 
parameter matrices) vary in magnitude with the frequency. How- 
ever, the diagonal elements of each of these matrices separately 
will be equal. Similarly, the off-diagonal elements of each 
of these matrices will also be separately equal. 


e.g. 


a p p 
pap 
-P p a 


For a 3x3 matrix of such a struct tire, there exists a 
frequency independent eigenvector matrix (vide Section 2.5) 


110 
10 1 
1 -1 -1 
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Hence for the case of an ideally transposed double 
circuit line» the matrices CQ^i] and [0^2^ equation 
(5.8) take the frequency independent form 


[Qil] = [0^2] = 



0 

1 

-1 


Thus a transformation matrix, which decouples the 
[P] matrix for the case of the transposed double circuit 
line at each frequency, inspite of the variation of [P] 
matrix with frequency, is given by 

1 1 0 - 1-1 0 

1 0 1-1 0-1 

1 -1 -1 -1 1 1 
110 110 
1 0 1 1 0 1 

1 - 1-1 1 -1 -1 

3.2.2 Untransposed line: 

For the case of an untransposed double circuit 
line# the eigenvalue, eigenvector analyses [21] for the 
matrices ([P^l+Cl^,]) and {[Pg^l-CP-i,]) are to be performed 
at each frequency to obtain the modal transformation 
matrix [Q]» 



3.3 CONCLUSIONS 

This chapter describes the development of a modal 
transformation matrix for a double circuit transmission , 
line of mirror symmetric conductor configuration with 
respect to a vertical plane. If the double circuit line 
is ideally transposed, it has been shown that a constant 
modal transformation matrix, independent of frequency, 
exists. This reduces the computational effort,' needed for 
the calculation of transients in the case of a transposed 
double circuit line. This model has been made use of for 
analysing energising and fault transients on double circuit 
lines in the subsequent chapters. 



CHAPTER 4 


EAULT IHITIaTIHG TRANSIENT ON A 
DOUBLE CIRCUIT TRANSMISSION LINS 

4.1 INTRODUCTION 

Switching overvoltages on an EHV transmission line 
arise due to energisation and re-energisation of the line. 
In practice, it is desirable to contain these overvoltages 
to a value as low as possible in order to reduce the cost 
of insulation. The various control methods such as the 
multistep resistance switching, controlling the timing of 
the closure of the individual poles of the cir.cuit breaker 
etc. are effective in restricting these overvoltages to a 
value as low as 1.5 p.u. When such effective methods of 
controlling overvoltages are employed, the overvoltages 
that arise due to fault initiation may become a limiting 
factor. Ninety percent or more of the faults, experienced 
by a power system are single line to ground faults. In 
the event of a single line to ground fault, transient over- 
voltages of the order of 2 p.u. may be experienced by the 
sound phases. This calls for an accurate evaluation of 
transients on a transmission line due to single line to 
ground fault initiation. 

The transient due to fault initiation on single 
circuit transmission line has been studied by various 
authors. Zimbark and Legate [22] have used the TNA to 
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evaluate the fault transient and have presented a lattice 
diagram approach for a theoretical study of the prohlem 
with certain approximations, later exhaustive analogue 
computer studies of the problem have been reported [23,24]. 
Boonyubol et al [ 25 ] have made use of the laplace transform 
method and have calculated the fault transient, employing 
the residue theorem for the inverse Laplace transformation. 

In all the above studies, the lines were assumed to be 
ideally transposed and the line parameters, computed at a 
single frequency have been assumed to remain constant in 
the whole frequency range. 

The double circuit transmission line is a common 
feature of any power system. In this chapter, the appli- 
cation of modified Pourier transform method for studying 
the transients due to fault initiation on a double circuit 
line is described. A typical double circuit line has been 
considered, taking into account the mutual coupling between 
the two circuits. A single line to ground fault is applied 
on phase ’a' at any point along its length and the resulting 
overvoltages on the sound phases are calculated. Both the 
cases of an ideally transposed line and a completely 
untransposed line are considered. For the case of an ideally 
transposed double circuit line, a frequency independent modal 
transformation exists, as shown in Chapter 3 and hence the 
calculation of the transient for this case involves much 
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less computation. The frequency dependence of line para- 
meters is included, using Carson's formulae [18] (Appendix B) . 

4.2 SIMUIATIOIJ OP PAUIT 

The principle of superposition has been applied for 
simulating fault on the transmission line. The system 
considered is linear and therefore superposition can be 
applied to the problem. Initially before the fault occurs, 
the line has been assumed to be under steady state. 

The steps involved in the simulation of fault can 
be explained with the help of Figure 4.1 as follows: 

(i) First assuming the fault not to be present, the 
voltage at the fault point and the required response at any 
point of interest on the line are determined as functions 
of time. As the line is in steady state, these functions 
are pure sinusoids at power frequency. 

(ii) As the second step; the voltage at the fault point 
on the faulted phase, calculated in step (i) is injected as 
cancellation voltage with reversed polarity, across the 
faulted phase and ground, at the same point. The source 
voltages at the ends of the line are short circuited. How 
the transient response of the line to this injected voltage 
alone is calculated as a function of time, starting from 
the instant of fault initiation (taken as zero time). 

(iii) The superposition of these two responses gives 
the total response as a function of time, subsequent to 
the initiation of the fault. 




flGf'RB 4.1.2 


PRIBCIFLE OF SltPERPOGI'TTCN VZBD FOR FAULT 


1,1 
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The steady state voltages to be evaluated in step (i) 
are nearly 1 p.u. The calculation procedure for the 
evaluation of the transient response to the injected voltage 
of step (ii) is described in the next section. 

4.3 evaluation OE TRANSIENT RESPONSE FOR THE 
INJECTED VOLTAGE 

The general solution of the transformed transmission 
line equations for the transformed phase voltage and current 
vectors at any point, x units away from the source end is 
given Ly 

V(x,w) 5= [Coship x] ^ 4 - [Sinh ^px] B (4.1) 

I(x,(o) = -[Z]”^[i)i]{ [Cosh ijix] B+ [Sinhr|,x] 1} (4.2) 

where [i|/] = [QILyIEQ]”^ 

[Cosh i|)x] = [Q][Cosh Y3t][Q]’"^ 

[Sinh ^|;x] = [Q][Sinh Yx][Q]"‘^ 

[y] = the diagonal matrix, whose diagonal elements 
are the square roots of the eigenvalues of 
the matrix [P] = [Z][Y] 

[q] = the eigenvector matrix (modal transformation 
- matrix), which maikes [Q] ^[P][Q] diagonal 

and A and B are arbitrary constant vectors, which have to 
be evaluated by forcing the terminal conditions of the line. 

Referring to Figure 4.1» let, the fault be located at 
a distance X units from the left end. The line section to 
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the left of the fault point F is considered first. The 
terminal conditions at the left end of the line are as 
follows : 

(i) As the two circuits are paralleled at this end, 
the transformed voltages of the connected phases of each 
of the circuits at this end, obtained by substituting 

X = Z in the general solution expression (4.1) should be 
equal. 

( ii ) The transformed phase cuarrents through the 
terminating source impedances at this end are the sum of 
the transformed currents fed by the corresponding phases 
of both the circuits, calculated by substituting x = X in 
the general expression for current (4.2)'. Further, the 
transformed voltages and currents at this end of the line 
should satisfy the transformed voltage-current relationship 
of the terminating impedances. 

These conditions can be put down mathematically in 
the following form: 

Let [CH] = [Gosh ij^Z] 

[SH] = [Sinh ipZ] 

[SC] = -[Z]“^ [4»][Cosh afX] 

[SS] = ~[Z]“^ [![;][ Sinh 

Then denoting the elements of the above matrices 
by lower case letters in double subscript notation 
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ohii 

°^12 



sh^ 

shi 2 

* • * ®^16 

°*"21 

Ohgg 

« • • oii^ ^ 

A + 

six^ ^ 

sli 22 


ohji 

ch ^2 

1 

• • • Cilry^ 

. i 

fhi 

sli ^2 

• • • sli^^ 


( ssj^]_+ss5j_) . 


(asig+ssgg) 


(ssjg+ss^g) 

(3336+83^6) 


A 


(sCii+SCgi) 
( 8021+305^) 

(sc^i+sc^l) 


. (sc^^g+scgg) 
. (sc 26 +sc 5 g) 


bL ( a+o w)l 


B 



““"ei 

0hg2 

... ., chgg 


shgi 

six^ 2 • • • . 

^‘‘ee 

= 


ohj2 

. . . ch^g 

M 

sh^i 

slx^2 • • • 

^*^56 


“‘‘41 

oh^2 

• * * *^^46 

' 

®^41 

slx^2 • 

31^46 


B 


(4.3) 


J 



^iXjU^L oh. 1 ^ • • • otx ^ 


^^11 ^^12 • * * ^^16 

= 

1 ^^2 2 • • • 6 

A+ 

8^2 ^ 2 ^ * ^^2 6 


oh.^^ oh^2 9 » ^ olx^^ 


sh^ ^ six^2 • • • 


B 


(4.4) 


where 1 is the source inductance, 
s 


Prom (4.3) and (4.4 )» B and A can be related by 


the expression, 
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1 = [Tj] A 


( 4 - 5 > 


wh. 6 r“e ~ 


[P 


xa^ 


[l 5 hii-(a+D(o)ig(sc^l+scg ] [ah^g-(a+d(D)I»g( sG;L 6 '^sc^g )] 

f[sh 2 i-(a +3 ti^lg(sc 2 i+sc^ ]-. •[ 56 

[sh^^-(a+jw)Ig(sc^j_-+sc^ 2 ^]^. *[ sh^g--(a+j'^ )LgC sc^gf-sc ^^)I 
( shg 3_-sh^^ ) - . . ( sb-gg “Sli^g ) 

(sli^^-sh 2 ]_) •> • (sh-^ 5 -sh 26 ) 

( Sk^g-Sh^g ) 




^ 3 ^ 1 = 


|[(a+j(o)I»g(ss^^+ssg3_)-otij_J .. .C(a+3o3)ig(ss^g-4ss^g> 

<a+3w)Ig(ss2i+ss5L)“Cti2 J .- £( a+D w) Is(ss26-+ss5g>-ch2g ] 

(a+ 3 w)lQ( ss 5 -L+ss^j_)-cln^ 2 ^ .-.E( a+j w) ^Lglss^g 4 ss^g>-cli 3 ] 
( ch^Li-chg 2_> - . C cti;L6“°'^66 ^ 

(ch2i-cli^l) •-* C°^26""°^56^ 

(cli^3_-cli^]_) • *• 


F*roiii tHe general solution espregsiorus ( 4 .l) aad ( 4 - - 2 ) axi. 
equia-ti*on ( 4 . 5 ), the input adiitifctsnce of the left aecthar 
of t»ie line in frequency do na-in, ¥hie n vie¥ed from tine fawlt 

poSri is, 

-Ir 


[yJ = -[Z] ['!'][%] 


( 4 . 
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37 ajiioopittia.i^ ; a a simi lar procedure, t:he input admittance 
oC ~th«3 ni -glut sjeiectitiMiiiiiof the line in the frequency domain, 
Tlewead fr‘‘oicQ tftb e : fa'aiiLJQt point can he proved to he 


!Z 


ir-l- 


j][r 


(4.7) 


wfcie Te ['l.,vTr] is-S' olxnteai.ned hy substituting { I - X) in the 

I Utma 1 

pdELcas of I iin- tline escipression for 

Ihae to ot naLDliitiHipiut admittance of the line aa viewed 
firom thie fsamlto pjoomnit [T] = [T^] + • 


looir-ajit thcQie fsraalt point, the transformed voltage on time 
raalfceci ;jh._asse,«, Wlj .('line modified Fourier transform of the 
Lnjectoei ir'oli.taEijeea aoioiirrce) is known and the injected currents 
o n time s* oiun® Hilis-Msses b ar e z er o . 


P'orrcifiiiSig 1 1 nes ae conditions, we get 


[-.ID 


- — 


— — 



I 

a 



0 

■'^0 

=: 

0 

^0- 

, 

0 

h’ 


0 



0 


(4.8) 


I&*oiiintdfclie&eBal)Ooowe equation, the transformed voltages 
on tliae sso- umd pliftaa.sesss at the fault point are given in terms 
oe «liBe Iliuorwi! jDuaijtbfcy Vp hy the expression 
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^22 ^23 ‘ • * ^26 

-1 . 

^^21 



^32 ^33 ^”36 


^31 


= 

^42 ^43 • * ' 5^46 


^41 

^b' 


7^2 ^53 • * * ^56 


^51 

V , 


^62 ^63 ^66 

1 

^61 


The transient response, given above in the frequency 
domain is then transformed hack to the time domain by perform- 
ing the numerical integration. 

4.4 TERIPICATIOJi OP THE GOMPOTm IROGRAM 

To check the computer program, developed for the 
fault study, the single circuit line examples, studied by 
Boonyubol et al [25] is solved. The data of the example 
is given in Appendix C. 

The voltage waveforms at the fault point on the 
sound phases are plotted in Figure 4.2, These results agree 
well with those of Boonyubol et al. 

4.5 SYSTEM USED FOR THE STUDY 

A typical 400 KV double circuit transmission line 
of configuration [lO], shown in Figure 4.5 is next analysed 
for the evaluation of fault initiating transient. The line 
is fed on both the ends by similar sources. In the example 
considered, both' the source sides are represented puredy by 

lumped impedances. The values of these source impedances are 

i. I. 1. i- / 

i CENTRAL IMRAm' 
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A\. o5 = 



54 


derived from the system short circuit power and voltages 
at these buses- This simplified model for the source is 
somewhat approximate- However, here the aim is just to 
illustrate the application of the modified Courier transform 
method of fault transient analysis for a double circuit line 
and hence source sides are not rigorously represented. For 
any particular system study, the source representation can 
be easily modified to suit particular requirements. For 
the present thesis, more complicated source representation 
is not considered to be necessary. However, there is no 
limitation for the applicability of the method for a 
detailed source side representation. The other lines 
terminating at the ends of the line under consideration 
can be represented quite accurately by equivalent shunt 
elements at these buses, the values of which are , obtained 
by finding out the input admittance offered by these lines 
at these buses, taking into consideration the distributed 
nature of these lines and the terminal conditions existing 
at the remote ends of these lines. 

The line parameters are calculated using the Carson's 
formulae. The effect of the presence of the ground wire is 
included by making the necessary correction to the line 
parameter matrices. The 50 Hz line parameters and other 
system data are given in Appendix D. 
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4 .6 COMPUTATIONAI. RESULTS 

4.6.1 Fault Overvoltages on Single Circuit ard 
Double Circuit Lines ■: 

The transient due to initiation of a single line 
to ground fault on the double circuit line, as described 
in the last section, is studied. In this study, the line 
is assumed to be ideally transposed and the 50 Hz line 
Parameters are used in the calculations. The fault is 
considered to be at the mid span on phase ' a* . It is 
assumed that the fault occurs at the instant when the 
voltage of phase 'a’ at fault point is at positive peak. 

Next, the presence of the second circuit is ignored and the 
fault study is conducted for the single circuit line, keeping 
all the parameters the same. A comparison of the waveforms 
of the sound phase voltages at fault point for both the cases 
of single and double circuit lines is made in Figure 4.4- 
A maximum peak of -1.78 p.u. is observed to occur on 
phase- c' of the second circuit for the case of the double 
circuit line representation. For the single circuit line 
case, the maximum peak voltage of -1.98 p.u. occurs on 
phase ' c' . The voltage waveforms for the double circuit 
line case clearly exhibit more high frequency oscillations 
than those for the single circuit line case. This appears 
to be due to the fact that there are three distinct modes of 
propagation with three distinct velocities of wave propagation 
in the case of transposed double circuit lines; whereas for 
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the case of a single circuit line, there are only two 
such distinct modes and velocities of propagation. Among 
these three inodes for the double circuit line, two are 
ground inodes, corresponding to the sum ( [Pg^]+[P^] ) and 
the difference propagation coefficient matrices. 

The other one corresponds'to the aerial mode of the sum and 
difference propagation coefficient matrices (vide 
Section 3*2). 

4.6.2 Transposed versus Untransposed Line: 

The fault study is next conducted for the double 
circuit line, now assuming the line to be untransposed. 

For this case, the modal transformation is evaluated at 
each frequency of the numerical integration range. The 
line parameters ore again assumed to be the 50 Hz values. 

The phase voltage waveforms at fault point (mid span) for 
the untransposed case are compared with those of the 
transposed case in Figure 4.5* It is to be noted that 
the maximum peak voltago in phase - c' is -1.97 p.u. for 
the untransposed configuration, whereas the maximum peak 
for the transposed configuration is -1.78 p.u. 

4.6.3 Effect of Incorporating Frequency Dependence 
of Line Parameters: 

Next the continuous variation of line parameters 
with frequeixjy is included in the calculation and the 
fault studies are repeated for both the transposed and 
the untransposed configurations of the line. The voltage 
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wavefonas for the cases of transposed and untransposed 
lines are plotted in Figures 4.6 and 4.7 respectiveHy . In 
each case the 'waveforms with and without incorporation of 
frequency dependence of line parameters are shown. The 
reduction in the maximum peak voltages when frequency 
dependence is included are about 11 percent and 13 percent 
for the transposed and untransposed cases respectively. 


The maximum peak voltages, obtained in the various 
cases for mid span fault on the double circuit line are 
summarized in Table 4.1. 


Table 4.1 


Case 


Max. peak 
voltage p.u. 


Transposed line, without frequency dependence -1.78 
Transposed line, with frequency dependence -1.585 

TJntr anspo sed line, without frequency dependence -1.97 
Untransposed line, with frequency dependence -1.72 


4.6.4 Overvoltages for Off-Centre Fault locations: 

So far in this chapter, the studies conducted 
have been for the mid span fault location. Next the case 
of the transposed line is considered and the study is 
repeated for different intennediate fault locations, ignoring 
the variation of line parameters with frequency. The sound 
phase voltage waveforms at fault point for the case of fault 












at a distance of 120 kms from one end of the line are plotted 
in Figure 4.8, The maximum peak voltages for various locations 
of fault points are given in Table 4.2. It is to be noted that 
the transients for fault locations symmetric with respect to 
the mid point of the line are identical because of the symme- 
trical nature of the system considered. 


Table 4.2 


Distance of the fault point 
from the line end in kms. 

(Total line length considered 
is 480 kms.) 

Maximum peak 
voltage in 

p.u. 

240 

-1.78 

180 

-1.57 

120 

-1.5$ 

60 

-1.40 

8 

-1.36 


The above off-centre fault location studies indicate 
that the maximum overvoltage occurs in the case of the mid 
span fault . 

In all the above studies, the fault is assumed to 
occur only on the bottom most conductor (Figure 4,3) at the 
instant when the voltage of that phase at that point fs passing 
through its positive peak. For the case of untransposed lines, 
if the fault occurs on one of the other two phase conductors, 
the fault transients and the associated maximum peaks will be 







different. However, there is no restriction on the 
applicability of the developed algorithm for studying the 
transients due to fault on the other phase conductors* As - 
the aim of the thesis is to present a method of analysis 
using the modified Fourier transform method for double 
circuit lines, rather than any particular system study, 
consideration of various cases of data permutation have 
not been studied. It may be pointed out here that the 
case of considering the fault on the bottom-most conductor 
need not be the most serious from the point of view of 
determining the system insulation requirements. Such 
considerations are outside the scope of this thesis, but 
are essential in carrying out actual system studies using 
the method proposed in this thesis or using any other 
method of analysis. This comment also applies for other 
examples that are considered in subsequent chapters. 

4.7 CONCLUSIONS 

The transient due to initiation of a single line 
to ground fault on a double circuit transmission line has 
been studied in this chapter. A mathematical analysis of 
the overvoltage transients, developed on the sound phases 
following the fault initiation at any intermediate point 
on the line has been presented. The double circuit line 
has been accurately modelled for the study, taking into 
account the mutual coupling between the two circuits. 

The study makes use of the modified Fourier transform method. 
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Both the cases of transposed and untransposed 
configurations of the line have been considered. As a 
frequency independent modal transformation exists for the 
case of the transposed double circuit line* as shewn in 
Chapter 3»- this case involves much less computational 
effort. In the case of the untransposed line, the evalua- 
tion of the modal transformation matrix at each frequency 
in the numerical integration range is required. The maximum 
peak voltage in the case of the untransposed line is observed 
to be nearly 10 percent higher than that of the transposed 
case. 

The effect of frequency dependence of line parameters 
on the fault transient overvoltages for both the transposed 
and untransposed lines is investigated , using Carson's formulae 
for the calculation of line parameters. The inclusion of 
frequency dependence of line parameters is observed to reduce 
the maximum peak voltages by nearly 11 percent and 13 percent 
for the cases of transposed and untransposed lines respectively. 

Finally the fault location is changed to various off- 
centre points along the line and the fault transients are 
studied. The mid span fault is found to give rise to the 
maximum peak overvoltage. 



CHAPTER 5 


EHBRGISING TRANSIEHT OH A DOUBLE CIRCUIT 
TRA.HSMISSION LIRE 


5.1 IHTRODUCTION 

An EHV transmission line is subjected to high, over- 
voltages because of energisation and re-energisation of the 
line. An accurate evaluation of these transients is 
essential for the economic design of the insulation. A 
double circuit line being a common element of any power 
system, a study of the transients due to its energisation 
is of importance. In the existing literature, the double 
circuit line transients have been analysed mainly using the 
transient network analyser [26]. In this chapter, the 
application of the modified Fourier transform method for 
analysing these transients and the results of these studies 
for both transposed and untransposed configurations of a 
typical line axe presented. 

5 .2 PROBLEM DESCRIPTION 

In practice, the energisation of a double circuit 
line is done in two stages. First one circuit is energised 
and after an elapse of time, the second circuit is closed. 
There is bound to be a delay of a few seconds at least if 
not a few minutes between the closing of the first and 
second circuits of a double circuit line. During this delay 



period, the transient arising in the first circuit due 
to its energisation would have died down and this circuit 
would be in steady state. Thus essentially two transients 
are involved in the study of the energisation phenomenon 
of a double circuit line, viz. (i) the transient due to 
the closing of the first circuit, while the second circuit 
is unenergised (ii) the transient due to the closing of 
the second circuit when the first circuit is already in 
steady state. 

5.3 FIRST CIROUIT ERIR&ISATIOIT 

In this section, a mathematical development of the 
analysis of transient due to the first circuit closure of 
a double circuit line, open at the receiving end is descri- 
bed » 

As described in the earlier chapters, the general 
expressions for the transformed voltage and current vectors 
at any intermediate point on a double circuit transmission 
line at a distance x units away from the sending end are 
given by the expressions 

V = [Oosh i|>x] E + [Sinh ipx] M 

1= { [Coshii> x] M+ [Sirih T{jx] K} 

where K and M are constant vectors, which are determined 
by the terminal conditions of the line and [\|>]» [z], 

[Cosh and [Sinh i|<x] are the same as defined in the 

earlier chapters. 
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If the voltage transforms at the sending and 
receiving ends of the line (V and are known, these 
conditions when forced oh the general expressions determine 
the constant vectors K and as follows: 

K = 7 

~ — s 

M = [Cosech it - [Coth i^J,] V 


where I is the length of the line. 


Hence the sending and receiving end transformed 
currents can he expressed in terms of the transformed 
voltages at the two ends as 


[A] [B] 


[B] [A] 


(5.1) 


where [a] = [ [Coth 

and [B] = -[z]”^[t{;][ C osech . 

1a] [b1 

In the above expression ^ ^ ^ can be inter- 

..[B] Ca]| 

preted as the short circuit admittance parameter matrix of 
the transmission line, treated as a two port network. 


How considering a double circuit line, which is 
connected to a voltage source with source resistance Rg 
and source inductance at the sending end and whose 

O 

receiving end is open (Figure 5.1*1), the voltage source 
with source impedance at the sending end can be replaced 
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by an equivalent source with shunt admittance as shown in 
Figure 5.1.2, 

As the second circuit of the double circuit line 
remains unenergised for the case under consideration, the 
forcing function is converted as a current source so that 
xij-hen the current injected to the second circuit phases at 
the sending end are talcen as zero, the problem is directly 
simulated. 


Now the double circuit line with the source admittance 
at the sending end can be viewed as a two port network and the 
incorporation of the shunt admittance at the sending end just 
introduces this additional admittance to the self admittance 
term of the sending end port of the short circuit admittance 
parameter matrix of the line. Therefore for this case, the' 
two- port admittance equation takes the form 


X' 


([Al+CTg]) 

■[B]’ 


V 

~s 

0 


[B] 

[A] 

1 


A 


where [a] and [B] are as defined in equation (5.1) except 
that they are now (6x6) matrices as we are dealing with the 
double circuit line. 


Jlere 



0 . 

0 

0 
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l/[Rg+(a+ja))l^] 0 0 0 0 0 

0 l/[Rg+(a+0a))lg] 0 0 0 0 

l/[ll 4-(a+ja))Ii ] 0 00 

o S 

0 0 0 0 

0 0 0 0 

0 0 0 0 


0 

0 

Equation (5.2) can be manipulated to write and 
V in terms of I as given below. 

4 = { [B]-([A]+[rg])[B]‘^[A]r^ Ig 

Is = ■ { [A]+[rs]-[B][A]-^[B] )-! 

The calculation of matrices [^], [Cosech 
[cothi(;jt] and hence the calculation of matrices [a] and [B] 
differ for the transposed and untrans posed configurations of. 
the double circuit line. It has been shown in Chapter 3 that 
because of the special structure of the propagation coefficient 


= 


0 

0 

0 

0 


0 

0 

0 

0 


and 


la = 


f Rg+( a+ jw ) Ig ] 
a+0 w)l^] 


s- 


n 
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matrix for the case of the transposed double circuit line, 
a frequency independent constant modal transformation 
matrix exists for this case. This very much simplifies the 
calculation procedure for the transposed configuration. 

For handling the untransposed double circuit line, the 
eigenvalue eigenvector analysis of propagation coefficient 
matrix is to be performed at each frequency to arrive at the 
modal transformation matrix at that frequency. The frequency 
dependence of the line parameters can be directly incorpor- 
ated in the analysis by using the appropriate line parameter 
values at each frequency, calculated using the Carson's 
formulae. 

5 .4 GOMFOTATIONAji RESULTS 

[Rie results of the study of single circuit energi- 
sation of a typical 400 E7 double circuit line, open at the 
receiving end are presented in this section. The configu- 
ration and parameters of the line are the same as for the 
example considered in Chapter 4 (Appendix D) for fault 
studies. The length of the line is 240 kms and the source 
is taken as purely inductive of value 0.1 Henry/ phase. The 
effect of the presence of the earth wire is included by 
making the necessary correction to the phase parameter 
matrices of the line. The earth resistivity is assumed to 
be 100 ohm-meter. 
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5.4.1 Transposed Line Versus Untransposed Line: 

Tlie phase voltage waveforms at the open receiving 
end of the first circuit for the transposed and untransposed 
cases are compared in Figure 5.2. For both the cases# the 
energisation has been done when phase *a’ source voltage is 
at positive peak. From the plots, it is observed that a 
maximum peak voltage of -2.45 p.u, occurs on phase ’a' at 

9.2 m.secs. for the case of transposed configuration. The 
maximum peak is 2.24 p.u. for the case of untransposed 
configuration and this occurs on phase 'c* at 12.8 m.secs. 
Because of the low order of magnitude of induced voltages 

on the unenergised second circuit, they have not been plotted. 

5.4.2 Single Circuit Versus Double Circuit Line 
Re pr e se nt a t io ns ; 

Next a set of energising studies are conducted, 
neglecting the presence of the second circuit with a view 
of comparing these results with those of double circuit 
representation. 

For the case of transposed line, there is no appreci- 
able difference in the waveshapes of receiving end voltages 
for the single circuit and double circuit representations. 

The reason fer this appears to be that the induced charges 
on the three phases of the unenergised second circuit 
cancel with one another even during the transient period 
in the case of the transposed double circuit line. 
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Th.e VO linage waveforms at the open receiving end for 
the untransposed single and double circuit line represen- 
tations are compared in Figure 5.3» From the plots, it can 
be observed that the majcimum peak voltage of 2,22 p.u. 
occurs on phase 'c' at 13.26 m.secs. for the case of single 
circuit representation; whereas a maximum peak of 2.24 p.u. 
occurs on phase 'c' at 12,8 m.secs. for the case of double 
circuit representation. 

5-5 SIMULATION' OF SECOND CIRCUIT BNER&ISATION 

In this section, the analysis of the transient due 
to the closing of the second circuit, when the first circuit 
is in steady state, is described. Again the same example of 
a double circuit line, open at the receiving end and 
connected to^ a voltage source with source impedance at the 
sending end is considered and an algorithm for the analysis 
of the second circuit energising transient, using the 
modified Fourier transform method is developed. 

The closing of the second circuit is shnulated by 
voltage cancellation technique and superposition. The 
steps involved can be described with the help of Figures 
5 •4.1 and 5.4.2 as follows: 

(i) First assuming that the second circuit continues 
to be unenergised, the response at any point of -interest 
on the line ( particular Oy the voltage at the open receiving 
end) a nd the voltages across the circuit breaker terminals 
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of all the three phases at the sending end of the second 
circuit are calculated. At this stage, as the first circuit 
is already in steady state, this step of the calculation is 
a steady state analysis. 

(ii) The voltages across the breaker poles, calculated 
in step (i) are injected as cancellation voltages, with 
reversed polarity, across the same poles. The sending 
end terminals of the first circuit are short circuited to 
ground. The response at the above mentioned point of 
interest on the line is again calculated, now for these 
injected cancellation voltages alone. 

(iii) Superimposing these two component responses, the 
total response is obtained. 

Therefore, the procedure adopted for calculating 
the second circuit energising transient of a double 
circuit line essentially boils down to a steady state 
analysis, keeping the second circuit unenergised and 
finding the transient response. of the line for a known set > 
of sending end voltage conditions. 

5.5.1 Calculation of Steady State Hesponse; 

The method of calculating the steady state response 
at any point of interest on the line and the voltages across 
the circuit breaker poles at the sending end of the second 
circuit (step i of last section) is very similar to the 
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analysis of the transient subsequent to the closing of the 
first circuit when the second circuit is unenergised (des- 
cribed in Section 5.3)» fhe only difference is that now 
the transmission line equations become full differential 
equations and the variables to be solved for are voltage 
phasors instead of instantaneous values of voltages because 
of the steady state nature of the analysis* It is to be 
noted that now all the variables involved are phasors at 
source frequency and no other frequency is involved in the 
calculations. 

By adopting a similar solution procedure as in 
Section 5*3» the expressions for the sending end and 
receiving end voltage phasors in terms of sending end 
current phasors can be derived as follows; 

I 3 = ([A]+[T^]-[B][A]'^[B])"^ 

=f [B]-([A]+[r^])[B]-l[A]r^ I3 

where [a] = [z]~^[^ ][Goth ] 

[B] = -[Z]"^['|i][Ooseoh *1] 

I = length of the Hne 
[Z] = CR]+j Wg[L] 

[1.] = [q][y]Cq]"^ 

[y] = diagonal matrix, whose diagonal elements 

are the square roots of the eigenvalues of 
the matrix product [z][Y] 
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[Y] = u^Lo] 

Ug * angular velocity corresponding to source 
frequency 

CQ] = eigenvector matrix of the matrix product [z][Y] 
[Cothi|»Jl] = [Q][Coth YJl][Q]~^ 

[CosechtJjJ,] = [Q][Cosech 



V(Rg+3a^I-g) 

0 

0 

0 

0 

0 



0 


0 l/(R3+3m3L^) 

0 0 


0 0 
0 0 


0 

0 

0 

0 

0 

0 


0 0 

0 0 

0 0 
0 0 
0 0 
0 0 


Rg = resistance of the source/phase 

Rg = inductance of the source/ phase 

[R] = resistance/ unit length matrix of the line 
[L] = inductance/ unit length matrix of the line 
[C] =3 capacitance/ unit length matrix of the line, 

( -0.5+30.866 )/(R +3 w r ) 
(-0.5-j0.866)/(Rg+3u)glg) 


and 


0 

0 

0 



If the double circuit line is transposed, the matrix 
[Q] takes a simple form as shown in Chapter 3* For the case 
of untransposed double circuit iLine, the eigenvalue eigen- 
vector analysis is to be performed for the matrix product 
[Z][Y] to obtain the matrices [y], [Q] and Cq]“^. 


When the above analysis is performed for the 
double circuit line example under study in this chapter, 
the steady state induced voltages on the unenergised second 
circuit at the sending and receiving ends are calculated 
to be as follows: 


Transposed case; 

At receiving end - Phase-a* 

Hiase-b' 

Hiase-c’ 

At sending end - Phase-a' 

Ihase-b’ 

Ihase-c’ 

Untransposed case; 

At receiving end - Hiase-a' 

Phase-b’ 
Phase-c’ 
PhasPra’ 
Phase-b ‘ 
Phase-c' 


( -0. 83x10”^ -jO. 11x10”^) p.u. 
(O.igxlO"”^ -30.72x10“^) p.u. 
(-0.97 x 10*”'^+30. 19x10"^) p.u. 
(-0.28x1D"^+30. 77x10“^) p.u, 
(0.68 xlO""^ -30.2 9x10“^) p.u, 
(-0.22 x 10“^+30.80 x 10~^) p.u. 

(0.016+30.68x10”^) p.u. 
(-0.025-30.025) p.u. 
(-0.082-3*0.079) p.u. 

( 0 ,a]l,l+30. 40x10”^ ) p.u. 

(-O.O29-3O.025) p.u. 
(-0.083-3*0.077) p.u. 


At sending end - 
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The steady state voltages on the first circuit at 
the sending and receiving ends are calculated to he 1.03 P.u. 
and 1.07 p.u. respectively. 

Prom the above analysis, it can be observed that the 
order of magnitude of the induced voltages on the second 
circuit is small when compared to the corresponding first 
circuit voltages. Hence it can be assumed that no voltages 
are induced on the second circuit conductors during this 
step of the calculation. In that case, the voltages across 
the sending end breaker poles of the second circuit can be 
taken as the steady state voltages at the sending end of the 
connected phases of the first circuit in step (i) of 
Section 5.5. However, if .more accuracy is required, the 
induced voltages on the second circuit can be taken into 
account idiile calculating the voltages across the breaker 
poles of the second circuit at the sending end for the 
case of the untransposed line. 

5.5.2 Calculation of Transient Response for the 
Injected 7oltages: 

The step (ii) of this calculation involves evaluation 
of transient response of the transmission line for known 
sending end voltage conditions (viz. the first circuit 
terminals connected to ground and the injected voltages 
connected to the second circuit) with the receiving end 
terminals kept open. 
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These conditions, when substituted in the general 
solution expressions for transformed voltage and current 
vectors at any point on the transmission line lead to the 
evaluation of the constant vectors K and ^ as given below: 

K = V 
— — s 

[CoshT/jJl] M + [Sinh\(;£] V = o 

«— g 

M = ■'[TanhifiJ,] 

Therefore the transformed receiving end voltage vector 
becomes 

m 

^ = CCoshr|;Jl] 7 - [Sinh [ Tanh\J;il ] 7 

= [ Sech 7 
[SechiJ»A] = [Q][Sech 

Cy] and [Q] are the eigenvalue and eigenvector 
matrices of the propagation coefficient matrix of the 
line respectively.. 

The inverse transformation of the transient response 
to the time domain is done by performing the numerical 
integration as described earlier in Section 2.2. 

5.6 COMPUTATIO MI RESULTS 

The same 400 K7 double circuit line example , 
considered earlier in Section 5.4 of this chapter, is 
taken up to study the transient due to second circuit 
closure and the results are discussed below. 


where 

and 



5.6.1 Transposed Versus TJntransposed Oonfigurations : 

The above study is done for both the cases of 
transposed and untransposed configurations of the line. 
Again the existence of the frequency independent modal 
transformation matrix for the transposed double circuit 

niuhss the computation for this case much easier when 
compared to the untransposed case. The comparison of the 
voltage waveforms at the open receiving end of the second 
circuit for the transposed and untransposed cases is done 
in Figure 5.5. For both the cases the closing of the 
circuit breakers is assumed to be at the instant when 
the source voltage of phase-a' is at positive peak. It 
is observed from the plots that a maximum peak voltage 
of 2,02 p.u. occurs on phase-a' at 1.02 m.secs. for the 
case of transposed line# whereas a maximum peak of 2.23 p.u 
occurs on phase-a' at 0.9 m.secs. for the case of untrans- 
posed line. The receiving end voltages on first circuit 
subsequent to the closure of the second circuit are not 
plotted because of the lower magnitudes of these voltages 
in comparison with those of the second circuit phases. 

5.6.2 Single Versus Double Circuit line Representations: 

Next a balanced three phase voltage of 1*03 p.u. 

(the steady state phase voltages at the sending end 
terminals of the first circuit before the closure of the 
second circuit) is applied to the second circuit and the 
resulting transient overvoltages at the open receiving end 
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are computed, neglecting the presence of the first circuit. 
This IS done to compare the transients for the single and 
double circuit representations of the line. 

For the case of transposed configuration of the line, 
there is no appreciable difference in the open receiving end 
voltage waveforms for single and double circuit representa- 
tions since the induced voltages on the second circuit due 
to the currents in the three phases of the first circuit 
balance one another. However, it is observed later in 
Section 7.7.1 that once the method is extended for handling 
the nonsimultaneous closing of the poles of the breakers, 
even for the transposed configuration of the line, the 
double circuit representation gives rise to a difference 
in voltage waveforms as compared to the single circuit 
representation. 

For the case of the untransposed line, a considerable 
difference in the waveshapes of the receiving end voltage 
waveforms is observed even for simultaneous closing of the 
breaker poles of the second circuit. A maximum peak over- 
voltage of magnitude -2,5 p.u. is developed on phase-a' at 
12 m.secs. for double circuit representation. On the*, other 
hand, a maximum peak of 2,02 p.u. occurs on phase-a* at 
1.26 m.secs. for the case of single Oircuit representation. 
Apar“t from the difference in peak values, the voltage waveforms 
(Figure 5.6) for the single and double circuit representations 
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differ widely in shape also. This seems to be due to the 
fact that for the case of untransposed double circuit lines, 
there are six distinct modes and velocities of propagation; 
whereas for untransposed single circuit lines, there are 
only three such distinct modes and velocities of propagation. 

5.7 EFFECT OF INOORIORATIOF OF FRSQTJENGI DEPEl'IDSNGE OF 
LIRE PARAMETERS 

Next the continuous vairiation of line parameter 
values with frequency is incorporated in the program and the 
voltage transients are computed for the various cases. How- 
ever, it has been observed that in all these cases, the 
effect of incorporation of frequency dependence of the line 
parameters on the voltage waveforms is not substantial. 

This appears to be due to the absence of the zero sequence 
ground currents in all these cases because of the simultane- 
ous closing of the circuit breaker poles and only -.zero 
sequence components of the resistance and inductance of the 
line vary considerably with frequency. It is observed 
later in Section 7.7.2 that when nonsimultaneous closing 
of the circuit breaker poles is incorporated in the program, 
the frequency dependence of line parameters affects the 
transient voltage waveforms considerably. 

5.8 CONCmiONB 

The modified Fourier transform method has been 
applied to analyse the energising transients of a double 
circuit line in this chapter. Algorithms have been 
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developed for computing the transient du^ to the closing 
of the first circuit with the second circuit unenergised 
and that caused by the closing of the second circuit when 
the first circuit is in steady- state. An example of a 
typical 400 KV double circuit line has been used to evaluate 
the transients for both transposed and untransposed confi- 
gurations of the line, IThe energising transients for the 
double circuit representation of the line have been compared 
with the transients caused by the energisation of the 
corresponding single circuit line (neglecting the presence 
of the other circuit). The effect of the inclusion of 
frequency dependence of line parameters on the energising 
transients is observed to be insignificant for the case of 
simultaneous closing of the circuit breaker poles. 



CHAPTER 6 


EXTSIJSIOH OP URMA Al© MILLER'S METHOD POR 
DOUBLE CIRCUIT LIl® EKERG-ISATIOH TRAHSIEHT 

ANALYSIS 


6.1 imODUCTIOI 

In this chapter, the theory of the Uram and Miller's 
method [ 4 , 5 ] of calculating the switching transients on 
transmission linos is extended for the case of double circuit 
linos. The existence of a constant modal transformation 
matrix for the case of transposed double circuit lines, as 
shown in Chapter 3 makes this possible. The example, tried 
in Chapter 5 (using the modified Pourier transform method) is 
taken up and the two transients, involved in the energisation 
of a double circuit line, viz., the transient due to the 
closing of the first circuit when the second circuit is 
unonergisod and that due to the closing of the second 
circuit, when the first circuit is in steady state are 
studied. Tho results obtained by both the methods are 
compared. As the Uram and Miller’s method is a time domain 
method, the frequency dependence of the line parameters 
cannot be directly incorporated in this method. The 
untransposed lines also cannot be handled by this method. 
However, it takes relatively less computer time. 
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6.2 EXTENSION OP IJRiiM AND MILLER'S METHOD FOR DOUBLE 
CIRCUIT LINE transients 

A brief outline of the Uram and Miller's method of 
solving electrical transients on transmission lines is 
given in Appendix E, In this section the theory is 
extended for the case of double circuit lines. 

The energisation of a double circuit line is to be 
studied in two stages. One transient is due to the closing 
of the first circuit t with the second circuit unenergised. 
The other possibility is that one circuit may be operating 
in steady state and the second circuit may be energised. 

Tho ovaluation of both these transients is essential to 
arrive at tho insulation level of the line and other 
equipments. 

It has been shown in Chapter 3 that a frequency 
independent constant modal transformation matrix exists for 
a transposed double circuit line. This makes the extension 
of tho theory of U!ram and Miller's method (described in 
Appendix B) possible for studying both the transients, 
involved in the case of a double circuit line. 

The equations (B.IO), (E.ll), (E.12) and (E.13) 
developed for the sending and receiving end voltage and 
current vectors in the Appendix E hold good for the double 
circuit line also. But now, the line parameters and hence 
[Z] and [I] become 6x6 matrices. The vectors k^(t), k 2 (t), 
2 t(t) and b(t) dre all of order 6. 
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The modal transformation matrix [T] takes the form 


[T] = 


C«] r[Q] 

!_£«] [Q]j 


where 



IS "the roohsl trs.nsfornifi'tl.oii msiiirix for 8. single cir 011111 
three phase line as described in Section 3.2.1. 

As described in Chapter 3# the line parameter matrices 
of the transposed double circuit line have the following 
special structures: 


II 


1%^ 

P 

[1] = 

1 — 1 







iy 



[C] 

ts 













In tho above matrices, each of the submatrices, [B. ] , 
[L ] aM [C ] is a 3x3 matrix, the diagonal elements of which 
are equal anti the off-diagonal elements of which are also 
equal. Bach of the submatrices [I»l 3 ] and [C^] is a 

3x3 matrix, the elements of which are all equal. 
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It has been shown in Chapter 3 that the double circuit 
line problem can bo reduced to two decoupled single circuit 
linos problem and the line parameter matrices of these two 
equivalent single circuits are ( [Rg^J+Lll^] ) , ( [l ]+[Lj^]), 
([C^]+[Cb]) and ([\]-[R^]), ( [ Cg^]-[G^] ) . 

• Denoting [R^]+[R^J by [R^], [l^] 4 -[l^] by [ij , 
CC^]+[C^] by [Cg], [R^J-ER^] by [R^], 
by [1^], [Cj-[0^] by [C^), 

the diagonal elements of [R ] , [l ] , [C 1 by R , L ,0 

s s s-* sd sd sd 

respectiveay; the off-diagonal elements of [Rg],[l ],[C J 

s s 

^^sod' ^sod' ^sod the diagonal elements of 

[R-(^] » [D^] » [C^] by ^dd' ^dd ^•^®P®otively and the 

off-diagonal elements of [R^],[l^],[C^] by R^^^, 

Cjjoa 3:‘ospoctivcIy, the matrix [ n] in equations (E.ll) and 
(E.I 3 ) (Appendix E) can be written for a transposed double 
circuit line as 


where 


Cn] * 



C^g] [0] 

lo] 

(W°a0>^ ° , 

0 0 
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0 0 



Co] = 


0 

0 



0 

0 

0 


0 

0 

0 


R^a = R ^ 

sO sd 

^dO “ ^dd 

A = R ^ 
30 3d 

^dO “ ^dd 

°a0 = °ed 

°d0 = °dd 


+ 2R 


sod 


+ 2P. 


dod 


+ 21 


sod 


+ 2X 


dod 


+ 20 


sod 


+ 20 


dod 


P 


P 


and 


^sl ^sd "* ^sod 
^dl = ^dd *■ ^dod 
\l " hd " ^sod 
^dl ^dd “ ^dod 
^sl = °sd “ “^sod 
^dl = ^dd ’ ^dod 


Hero it can be observed that R.t=Rt ;1,t=L-, 

dl si dl si 

*^dl “ ^sl showing that there are only 3 distinct" modes of 
pro pngn tion in the case of transposed double circuit lines 

6.3 FIRST CIROITIT .SHERGISITIOH TRANSIENT 

The first step in the Uram and Miller’s method is 
to express the vectors terms of the delayed 

functions a and ^ by forcing the boundary conditions, 
existing at the ends of the line. 
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Rcfarrins to Figure 6.1, which is a single line 
diegrh-^ for the first circuit energisation study of a double 
circuit line, open at the receiving end, the expression for 
the sending end current [Equation (E.ll)] can be written as 


-1 


[Q] 

-[Q] 


[«g]'^ [0] 


(fel)l “(^2^1 

i2 


[Q] 

[Q] 


[0] 


**^^2 ^2 


( 6 . 1 ) 


wh.Gr*o nnd ko tii'o pfirt/il/ioiiGd 3,3 


^- 1^1 


(^)l 


ard = 



“2 

1_ .--J 


The expression for the sending end voltage 
[equation (B.IO)] can be written as 




[Q] 

-[Q] 


^-1^1 ^-2^1 



[Q] 

CQ] 


(^ 1)2 + (^2^2 


As the second circuit is unenergised, forcing the 
coiidition 1^2 « 2 equation (6.1) yields 

i.e. (iCjL)2 “ (^2^2 ^'^(^2^1 ~ ^-i^l^ (6.3) 
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The other sending end terminal condition is determined 
by the soarco impedance. This can be written as 


CQ] { (ki)i + (^2)1 “ (^1)2 - (^ 2^2 ^ 


Expressing (|c^)2 in terms of as per equation (6,3), 

the above equation can be manipulated to the form 

1-1 ™ . 1 


111 ) 1 1 ^ 1 n 

- = [nd3+f"s3>](ki)a 




( 6 . 4 ) 


Ml _ 

g 


S 


S 


By forcing the condition that the line is open at 
the receiving end, wo get the relationship. 


k. “ n 


( 6 . 5 ) 


The deljiy functions ^ and k2 are related to the functions 
^1 and ^ as follows: 


c 


- ■'so^h 


rv ft - T )1 

^2 * 0 • L^12^^ ^sl'^-'l 





5 


^- 13 *^ - 1 i)]i 


= ■'*°dc/^do>‘ ,., 

'4 - ^.n)] 


iO'^2 


- o 


■(Rgl/2) •'(Cjl/lgl)* 


- t„,)] 


sl '-‘2 


-(R^i/a) /( 0 ^,/Vt 

® [il3(t - t^i)]2 


( 6 . 6 ) 


..•here T^o “Ityd-go ° 80 ) ■ "sl=»«l 3 i C 31 ) and 

■‘do °ao^' 

-(R^o/2) AO^q/I^o)* 


( IC 2 2^ ) ®®' o 


(k ) - 

'■‘^22^1 - ° 


- ^so) 
b2{t - 


-(H,l/2) /(O^/I ,)t 
(kgj)^^ = e bj(t - 


(6.7) 


( *1 )p 

«N»«' mmm 


(> 1 : 22)2 


-(Rdc/2>'''<WV‘ ,, 

° '>4(t - 


(k23>2 = ° 


.(R^j/2) /(Ogi/lgi)lt 

-(K 31 / 2 ) /(C3l/l3l> ‘ 


’’s*'*' ■ ■'si' 


bglt - T^^) 


¥ith tho help of eqoations (6.3) - (6.-7) an 
incromental solution for the four vectors k^, a and b 



at each instant of time and hence the response at any point 
of interest on the transmission line can be evaluated. As 
equation (6.4) is a differential equation, a numerical 
solution of tiiis equation becomes necessary to arrive at 
the values of at each instant. 

6.4 SSC0I3D CIRCITII} BN3RGISATI0N 

An algorithm for studying the transient due to the 
energisation of the second circuit, when the first circuit 
is m steady state is developed in this Section, The voltage 
cancellation technique and the principle of superposition are 
used for simulating the circuit breaker closing operation. 

The principle of superposition is applicable for this part 
of the analysis as the system considered is linear. 

The steps involved can be described with the help 
of Pigijres 6.2.1 and 6.2.2 as follows: 

(i) First assuming that the second circuit continues 
to be unonergised, the response at any point of interest 

on the line (particularly the voltage at the open receiving 
end) arxi the voltages across the circuit breaker terminals 
at the sending end of the second circuit are calculated. 

As the first circuit is in steady state, this step of the 
calculation is a steady state amlysis* 

(ii) The voltages across the breaker poles, calculated 
in step (i) are injected across the poles with reversed 
polarity as cancellation voltages and the sending end 



tcrainnls of too first circuit ore grounded. Now with 
thfso ion voltuges alone present, the response at 

the 'ibove fflentionod point of interest is again calculated. 

(iii) Swperiaposing these two conponent responses, the 

total rc;?p0:-i3c is obtained. 

The algorithm for finding the response to the 
injected voltage (stop ii) by Uram and Miller's method is 
developed as follows*. 

Referring to Figure 6.2.2, the terminal conditions 
at the sending end arc that the voltages of the first 
circuit arc aero and those of the second circuit are v^^^. 


Forcing these conditions and using the same 
notation, na used for the first circuit energisation case, 
equation (6.2) becomes 
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for i„ oquation (6.9) by equation (6.10), 

¥ € (yJ "t* 

(il)? =|tQr^Vi-{S2)3 (6.11) 

SubG’ti bii 1 1 Tig ( 6 # 11 ) in ( 6 • 10 ) , vo gsi 

li - (^)i (6.12) 

A3 the receiving end is open, from equation (E.I3), we get 

li = a (6.13) 

fhf*' functions k have been expressed 

Ir. terjua of delnjed functions and known quantities [Q]» 
by equations (6.11) - ( 6 .I 3 ) respectively. 

The delay functions a, and ^ are related to the 
functions and ^ by the sets of equations (6.6) and 
(6.7). Hence with the help of equations (6.6), (6.7)» 
(6.11), (6.12) and (6.13), an incremental solution for 
the desired response for the injected voltage can be 
c-alculatcd. 

6.5 COMPOTAflOMAli R«SUia?S 

The same example, tried in CSaapter 5, using the 
modified Fourier transform method is taken up and the 
Virm and Miller's method is used for evaluating both 
the first and the second circuit eiffirgising transients 
on the double circuit line. The modal components of the 
parameters for this line are given below. 



J-Oj) . 

O/raile ; = 0.00686016 H/mile ; 

0 = 0,954644:il:0“® f/mllo 

’'si “ ’‘'il “ 0-030330 n/milo ; a = O.OOM2587 H/mile 
O.J . C,J^ .-. 2.085753x10-® f/mile ; = 0.03 0 /mile 

I-dO = .00167236 li/iaile ; C^q = 1.620324x10“® f/nile. 

Tho ph'isc voltn^G waveforms at the receiving end of 
the flFBt circuit, subsequent to the closing of the first 
circuit :»ro plotted in Pigure 6.3. The breakers of all the 
three phria<j3 aaaumod to close simultaneously at the 
inr.tant when the aourco voltage of phase 'a' is at positive 
pe.-ik. A maximum pciak of 2.32 p.u. is observed to occur on 
phaoa *a' at 9.2 m.socs. for this case. The corresponding 
voltage WMVofoms, obtained using the modified Fourier trans- 
form method :\Tq also plotted in the same figure. It can be 
ohserv /d that there is a close matching between the waveforms 
obtairied by bOiih the methods. 

«RXt for the same line, the transient due to the 
closing oi’ the second circuit, when the first circuit is in 
steady jit.-to (described in Section 6.4) is studied.. The 
steady state 3 imilysis of Section 5.5.2 holds good here also 
and hence in step (ii) of Section 6.4» it can be assumed 
that no voltages are induced on the second circuit conductors. 
Hence, the voltages across the breaker poles of the second 
circuit (which is the injected voltage source) can be taken 
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as the steady state voltages a.t the sending end of the 
connected phases of the first circuit, Again the breaker 
poles of all the three phases of the second circuit are 
assumed to close • simultaneously at the instant, when the 
phase-a’ voltage is at positive peak. The phase voltage 
waveforms at the receiving end of the second circuit, 
subsequent to the closing of the second circuit are 
plotted in Figure 6.4. A maximum peak of 2,02 p.u. is 
observed to occur for this case at 0.9 m.secs. on phase-a' . 
These waveforms are almost coinciding with the corresponding 
ones, obtained using the modified Fourier transform method 
in Chapter 5. 

6.6 CONCIUSIONS 

The extension of the theory of Uram and Miller's 
method for solving the energisation transients on a trans- 
posed double circuit line has been presented in this chapter 
The same line, tried in Chapter 5 using the modified Fourier 
transform method, has been taken up and the transient due to 
the energisation of the first circuit alone, with the second 
circuit remaining unenergised and the transient due to the 
energisation of the second circuit, when the first circuit 
is operating in steady state, have been studied. The voltage 
waveforms at the open receiving end of the line, obtained 
using both the Uram and Miller's and the modified Fourier 
transform methods have been compared. The results by both 
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the methods have been found to be closely matching. As 
the Uram and Miller’s method is a time domain method, the 
frequency dependence of line parameters cannot be directly 
incorporated in this method. The untransposed lines also 
cannot be handled by this method. However, this method 
takes relatively less computor time. 



CHAPTER 7 


EHERCrlSATIOH TRAl^TSIEHT DUE TO SEQUENTIAL CLOSING OF 

CIRCUIT BREAKER POLES 


7.1 INTRODUCTION 

In reality, during tlie energisation of a transmissio: 
line, the poles of the circuit breaker do not close simulta: 
ously. However, such an assumption is generally made only 
simplify the analysis and to get an idea of the magnitude c 
the switching overvoltages involved. The difference in ti: 
in the closing of different poles of a circuit breaker is 
small, but this time is comparable to the time period of tt 
electrical transients in the overvoltages that develop, Th 
has been seen in earlier chapters to be of the order of a f 
m.secs. Therefore, for an accurate analysis, the non- 
simultaneous or sequential closing of the circuit breaker 
poles has to be incorporated in the program. The sequentia! 
closing of the breaker poles is bound to give rise to a hig 
order of overvoltages. There are basically two approaches I 
the existing literature for the incorporation of sequential 
closing feature in the Fourier transform method of studying 
the energisation transient on a transmission line. ¥edepohll 
and Mohamed [27] have proposed a method of solving the prob' 

by numerically evaluating the Fourier transform of the indt; 

i 

voltages on the unenergised phases, known at discrete ; 

intervals of time. On the other hand, the method proposed | 
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by Battison et al [28] solves the problem by using a 
convolution integral approach for taking into account the 
sequential pole closure. In this chapter, the method 
proposed by Wedepohl' and Mohamed is extended for the case 
of a double circuit line. After presenting the relevant 
theory for this extension, the double circuit line example, 
with receiving end open, of Chapters 5 and 6 is considered. 
The transient, arising due to the sequential closure of the 
circuit breaker poles of the . second circuit, when the first 
circuit is operating in steady state, is evaluated. The 
transient voltage waveforms obtained by the double circuit 
modelling is compared with those obtained by ignoring the 
presence of the first circuit. Also the effect of the 
continuous variation of line parameter values with frequency 
on the energisation transient for both transposed and un- 
transposed lines is investigated, 

7.2 INCORPORATION OR SEQUENTIAL CLOSINC OR CIRCUIT 
BREAKER POLES IN ROURIER TRANSRORM METHOD 

The sequential closing of the circuit breaker poles 
involves changes in the terminal conditions at the sending 
end of a t;ransmiseion line at the closing instant of each 
breaker pole. As the Rourier transform method is a frequency 
domain method and because the change in system conditions 
due to sequential closing occurs at different instants of 
time, an indirect solution of the problem becomes essential. 
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To start with, when only the f irst phase of the line 
is getting closed, the conditions to he forced at the sending 
end are that the currents on the other two unenergised phases 
are zero. Further, the voltage and ciirrent on the first 
phase at the sending end are determined by the source voltage 
and the impedance, present on this phase at this end. After 
a short time of the order of a few m.secs., the second pole 
gets closed, leaving only the third phase unenergised. At 
this stage, the current through the third phase is zero and 
the voltages and currents in the first and second phases at 
the sending end are determined by the source voltages and 
impedances present on these phases at the sending end. After 
the elapse of a few more m.secs,, the third pole gets closed 
and from this instant onwards, tbe voltages and currents on 
all the three phases are determined by the source voltages 
and impedances at the sending end. During the interval 
between the first and second pole closure, voltages are 
induced on the floating second and third phase conductors. 

In the Fourier transform method, these induced voltages at 
the sending end should be first calculated as functions of 
time and then the modified Fourier transform of these 
functions so- calculated at discrete interval of time should 
be numerically evaluated to be fed as input in frequency 
domain. Similarly during the interval between the second 
and third pole closure, voltage is induced on the floating 
third phase conductor because of the mutual effect between 
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this phase and the other two energised phases. The modified 
Fourier transform of this function in this tme interval, 
augmented hy the induced voltage on this phase during the 
interval between the first and second pole closure is to he 
numerically computed and fed as input voltage in frequency 
domain to this phase. 

The steps involved in the solution of the sequential 

closing transient can he summarized as follows”. 

(i) The first pole closes (say at zero time) and th.^ 
sending end boundary conditions at this stage are a known 
voltage on the first phase and zero currents on the other 
two phase conductors. Using the modified Fourier transform 
method, the system response at the point of interest 
computed as a function of time. In particular, the induced 
voltages at the sending end on the two unenergised phase 

conductors are also computed. 

(li) After a time interval q, the second pole closes 

and the new boundary conditions at the sending end are known 
voltages on the first two phases and aero current on the 
third phase. The transformed voltages to be fed as input 
on the second phase Includes the modified Courier transform 
of the induced voltage on this conductor in the interval 
0 < time < q (evaluated numerically) m addition to the 
transform of the source voltage, connected to this phase 
for time > (for which an analytical expression can be 

obtained). For these terminal conditions, the response 
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at the point of interest can be computed as a fimction 
of time. In particular, the induced voltage on the floating 
third phase conductor is also to be computed. These res- 
ponses are valid right from zero time. But the response, 
so obtained -will be a mere repetition of the corresponding 
response, obtained in the first step in the interval 
0 < t < T^. However, if it is computed again here, it 
can be utilized to get an idea of the error, introduced 
by the numerical integration and the numerical evaluation 
of the transform of the induced voltages involved in this 
st ep . 

(iii) After an elapse of time, finally the third pole 
closes at time T^. Now the transformed voltages on all the 
three phases at the sending end are known and hence the 
response at the point of interest for time > T 2 can be ; 

computed. Here the transformed voltages, to be fed at the 

sending end as input, includes in addition to the analytically: 
known transformed voltages on all the three phases, the 
numerically evaluated transforms of the induced voltage on | 
second phase in the interval 0 < time < T^ and that of the | 

induced voltage on the third phase in the interval I 

0 < time < T^, Again the computing of the response right 
from zero instant helps in giving an idea of the error, | 

i 

introduced by the numerical processes involved. 
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1,3 DEVELOPMENT OP ALGORITHM FOR HANDLING SEQUENTIAL 
CLOSING 


As described in the last section, the incorporation 
of sequential closing essentially involves finding the 
transient response of the line for three different sending 
end terminal conditions. The change in terminal conditions 
occurs at the closing instant of each pole. The terminal 
conditions are zero currents on unenergised phases and 
toiown voltages on the energised phases behind source 
reactances, which can be converted to equivalent current 
sources with appropriate shunt admittances. Hence, onoe 
a set of equations are formed for computing the transient 
response of the line for known sending end currents, it 
can be used repeatedly thrice with the suitable currents 
injected at each of the pole closing instants, to get the 
complete response. 


It has been derived in Section 5,3 that when a 
transmission line is viewed as a two port network, the 
transformed currents and voltages at the terminals are 
related by the expression. 


"is 


I — 1 

[B] 



}r 


1 — 1 t 

[A] 




where = [Z ]'”^£^3[Coth^il] 

[B] = -[Z]“^[i(»][Cosecht|iA] 

[Z] [R]+(a+ja))[L] 
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[Y] = (a+jw)[G] 

M = [QJEyILQ]”^ 

[Cothtf^ii] = [Q][Ooth 
[Coseoli4»Jl] = [Q][Cosecli y 

and [y] is the diagonal matrix, whose diagonal elements 
are square roots of the eigen values of the 
matrix product [Z][Y], 

When the open receiving end condition is forced into 
equation (7.1), and can be expressed in terms of 
Ig as follows: 

Is = ([A] - (7.2) 

Ir = ([B] - [A][B]"^[A])"^ Ig (7.3) 

When there is a shunt admittance [Y^] at the sending 
end, equations (7.2) and (7.3) take the form 

Is = { [A]+[^-[B][A]-1[B]}-1 (7.4) 

Ih ={[B]-([A]+[T^])[B]-kA]}~h^ (7.5) 

where now I_ includes the current drawn by the shunt 
admittance in addition to the current drawn by the line. 

In the equations (7.4) and (7.5), the quantities 

[^g] and Ig vary according to whether the different phase 
poles of . the breaker at the sending end is open or closed. 
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Initially, in the interval 0<t <l2_» when only the 
first pole is closed, [Y„] and are given hy the expressions 

S ““S 

l7TRs+(a+0w)lg] 0 0“ 

0 0 0 

0-00 


th] = 



0 

0 


where is the modified Fourier transform of the first 

phase source voltage and R„,l_ are the source resistance 

s s 

and inductance per phase respectively. 


In the interval < t < T2» first and second 

poles are closed and now 




l/[Rg+(a+ow)Ig] 0 0 

0 l/[Rg+(a+Ga))lg 0 

0 0 0 



Eg;j_/[^^g+(a+3 

0 


where E « is the modified Fourier transform of the source 

sa 

voltage of this phase (which becomes the line side voltage 
also on this phase from time onwards) [i.e. the modified 


Fourier transform of ^ 


and is the 
sa 
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numerically evaluated modified Fourier transform of the 
induced voltage on the second phase conductor at the sending 
end in the interval 0 < t < 

From the instant onwards, all the three poles are 
closed and hence 



|l/[Rg+(a+>)Lg] 0 0 

0 l/[Rg+(a+oa))lg] 0 

0 0 l/[Rg_+(a+ja))Lg} 

(Eg^+Vg^)/[l{g+(a+j w)Lg] I 


where is the modified Fourier transform of the source 

voltage of this phase (which becomes the line side voltage 
on this phase from time T2 onwards) [i.e, the modified 
Fourier transform of eg^C'fc ~ ^3^ ^ ~ ^2^^ 

the numerically evaluated modified Fourier transform of the 
induced voltage on the third phase conductor at the sending 
end in the interval 0 <t <12* 

The appropriate values of [T„] and when used 
in expressions (7.4) and (7.5) in the appropriate intervals 
of time, facilitates calculation of the desired transient 
overvoltage response in frequency domain and it can be 
converted to time domain by computing numerically the 
inverse transformation integral. 
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Here it can be observed that the numerical integra- 
tion of the transformed response to time domain can be limited 
to the extent of calculating the induced voltages at the 
sending end on the second and third phase floating conductors 
in the interval 0 <t<T^, subsequent to the first pole closing 
and calculating the induced voltage at the sending end on the 
third phase conductor alone in the interval T^< t < T 2 , 
subsequent to second pole closing. The transforming back 
of the receiving end voltages to time domain may be done just 
once at the end. 

The method of calculating numerically the modified 
Fourier transform of the induced voltages on unenergised 
conductors, known at discrete interval of time is given in 
the following section. 

A flow chart of the program for computing the 
energisation transient on a single circuit tranmission 
line with sequential closing of circuit breaker poles by 
modified Fourier transform method is presented in 
Appendix F. 

7.4 mJMERICAl EVALUATION OF THE MODIFIED FOURIER TRANSFORM 

For numerically evaluating the modified Fourier 
transform of a known function f(t), it can be divided into 
N strips of equal width, the interval of the strip being 
sufficiently small, such that the variation of the fimction 
within any of these intervals can be approximated by a 
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straight line. Then the modified Fourier transfoim of the 

function f(t) can be obtained by summing up the modified 

Fourier transforms of all these individual strips, A 

typical strip of the function and how it can be split up 

into components of step and ramp functions are shown in 

Figures 7,1,1 and 7,1.2. This strip can be split into 

two component step functions (one of magnitude 

starting at instant (n-l) At, the other of magnitude 

-f , starting after a delay A t at the instant n At) and 

A f^ 

two ramp components (one of slope , starting at 


A f . 


instant (n~l) At and the other of slope - , starting 

after a delay At at the instant nAt), 


When the modified Fourier transforms of the two 
component steps of all such strips are summed up, the 
transform of all such steps in the entire range of time 
is obtained as 




[f^ - f„ 

(a+jm) ' o N 


(7.6) 


The sum of the modified Fourier transforms of all 
component ramps can be given by 


F (w) = - 

At(a+jaj)^ 


N 


I Af [e”(^+3w)(n-l)At _ 


n=l 


g-(a+ja))nAt^ 

(7.7) 


An economy in the computation of F^(‘^) can be 
effected by expanding and grouping the consecutive terms 




PIG-TIRE 7.1j NEIffiRIGAI. EVAIiUAIIOH OF POTJHUR IRMSPORM 
OP A PUl'ICIIOR. 
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of (7.7) together. By this manipulation, equation (7.7) 
can he reduced to the form 

P ( m) = - — — ff —f —(f _f ) g“(a+ja))NAt 

^ At(a4-aa))^ ^ 1 ° ^ 

Hence the modified Fourier transform of the function 
f(t) in the whole range of time can be written as 


F(a)) = 


[f {1 - — 1 


(a+jw) ° At(a+j(o) At(a+ja)) 


+ { Illl f^(l + 1 __) } e 


-(a+3o3)NAt 


At(a+i) ) 


At(a+o m) 


H-l 


A t ( a +3 w) n=l 


r (f - 2f + -F p~(9-+0“)iiAt-| 

i-. ^^n+1 '^^^n + ^n-1'' ® -I 


7.5 SINGLE CIRCUIT LINE EXAMPLES 

7.5.1 Verification of Computer Programs 

To check the computer program, developed using the 
modified Fourier transfoim method, incorporating the 
sequential closing of circuit breaker poles, the transient 
voltage at the receiving end of a typical 400 KV single 
circuit line, computed by this method is compared with the 
results, obtained for the same example by Uram and Miller' s 
method. 
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This ■transmission line of Uttar Pradesh 


Stat, 


Electricity Boar'(3. is 400 kms, long and connects IatcIoiiow 
with Ohra. It h.as got two shimt reactors, one eaoh. at 
sending and receiving ends. The energisation stadv * 
done for energisation at Obra end, Figure 7,2 giv^Q 
system set up, system data are given in ^ 


The phase sequence of closing of the cireuii breaker 
poles is assumed for the study to be RBY and the swi-tching 
instants of B pha.se and Y phase are taken as 3,3 m 

♦ o S • 

and 6,7 m.secs, respectively (taking the instant: of closing 
of R phase as zero time). This set of switching instants 
are such that each pole gets closed when the source voltage 
of that phase is at peak value. 


The phase voltage waveforms at the receivijj^g 
obtained using both the methods are compared in Pig^^g Y 3 
There is a close matching between the waveforms, obtained 
by the two methods, A maximum peak of -3,01 p,u -i 4 .u 
case of the modified Fourier transform method and ..-2 88 p u 
in the case of the Uram and Miller’s method are observed to 
occur, in both the cases on phase R at 9 m,secs. 


This particular example contains reactors at 
both the ends of line additionally as compared to 
the system, considered in Section 7.3, while '^®'^elop^• 
the necessary theory for incorporating sequential closing 
feature in the modified Fourier transform method, Uhen 
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these reactors are introduced, the equations of Section 7.3 
get modified as given in the following para. 

in equation (7.1) can he 

interpreted as the short circuit admittance parameters of 
the transmission line. The introduction of the shunt 
admittances at both the ports just adds these 

admittance values to the self admittance terms of the 
short circuit admittance parameter matrix, for obtaining 
the overall short circuit parameters of the system. 

Hence equations (7.4) and (7,5), with the reactors 
introduced, read as follows: 

Vg ={ [A]+[Tgj^]+[Yg]-[B]([A]+[T^j^]3'[B]}“^ I 3 (7.8) 

(7.9) 


The matrix 


[A] [B] 

[B] [A] 


where 


l/(a'+D w)I'gj^ 0 0 

0 l/(a+j(D)Igj^ 0 

0 0 l/(a+ju))lgj^ 


and 


is the inductance of the shunt reactor. 


7.5.2 Effect of Incorporating Frequency Dependence of 
Line Parameters: 

In order to investigate the effect of frequency 
dependence of line parameters on the transient overvoltage, 
produced by the sequential closing of circuit brea.ker poles. 
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a single circuit of the double circuit line example, 
studied in Chapters 4,5 and 6 is considered,. This line 
is 150 miles long and has a source inductance of 0,1 H, 

The continuous variation of the line parameter values with 
frequency are calculated, using the Carson's formulae and 
are used in the calculations. The switching instants are 
assiimed to be 0, 3.3 and 6,7 m.secs. The phase sequence, of 
pole closing is assumed to be RBY, so that each pole gets 
closed when the source voltage of that phase is at peak 
value. The comparison of the transient voltage waveforms 
at the open receiving end with and without frequency 
dependence of line parameters is done in Figure 7.4. From 
the waveforms it can be observed that the maximum peak 
overvoltage value gets reduced considerably and the 
oscillations in the voltage waveforms get smoothened, when 
the frequency dependence of line parameters is considered, 

1 maximum peak of 2,72 p.u. occurs on phase Y at 9.3 m.secs. 
for the case without frequency dependence of line parameters; 
whereas a maximum peak of -2.36 p.u. occiirs on phase B at 
5.9 m.secs, for the case with frequency dependence of line 
parameters, included. The line is assumed to be uniformly 
transposed for both the cases, 

7.6 EFER&ISATION TEAHSIMT OH A DOUBLE CIROUIT LINE 

In this section, the necessary theory for studying 
the transient on a double circuit line, due to the sequential 
closing of the circuit breaker poles of the second circuit. 
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■when the first circuit is operating in steady state, is 
developed. Section 5.5 describes the method of calculating 
the transient due to the closing of the second circuit but 
for the simultaneous closing of the breaker poles only, for 
the energisation study, introducing the sequential closing 
feature also, the principle of superposition can be used, 
as shown in Section 5.5. But additionally here, the effect 
of induced voltages on the floating conductors prior to 
tbeiv energisation has to be considered. As described in 
Section 5.5, the procedure for calculating the second circuit 
energisation transient of a double circuit line involves 

(i) a steady state analysis, assuming the second circuit 
to remain unenergised and 

(ii) the determination of the transient response of the 
line for a known set of sending end voltage conditions. 

It has been shown in Section 5.5.1 that for the case 
of a transposed double circuit line, the induced voltages 
on the un energised conductors, obtained in the above mentioned 
steady state analysis (step i) will be negligible. Hence 
the voltages to be injected at the sending end on the second 
circuit conductors, for which the transient response has to 
be calculated (in step ii) are the steady state values of 
voltages of the connected phases at the sending end of the 
first circuit. The voltages on the first circuit conductors 
at the sending end for the calculation of the above transient 
response are zero. 
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The procedure for the calculation of this transient 
response at the point of interest (in this case the voltage 
at the receiving end) on the double circuit line may be 
described by referring to Figure 7.5. 


The general expressions for transformed voltage and 
current vectors at any point on a transmission line, subse- 
quent to a switching operation, using the same notation as 
used early in this chapter, are given by 

Y = [Goshtjixlg + [SinhifcclM (7.10) 

I = J { [Cosh + [Sinh (7.11) 

where K and M are constants, determined by the boundary 
conditions, existing at the ends of the line. 

In the interval 0 < t < T^^, only the first phase of 

the second circuit is closed. Hence in this interval, the 
voltage on the first phase alone gets the value of the 
injected voltage and the other conditions, to be forced 
at the sending end are zero ctirrents on the two unenergised 
conductors of the second circuit, in addition to the zero 
voltage conditions on the first circuit conductors. 


By forcing the known voltage conditions at the 
sending end, the first four components of K in equation 
(7.10) are determined as follows: 


== 0 


% - 
E. = B 


^2 = 0 ; 


K. 


0 


( 7 . 12 ) 


si' 


the modified Fourier transform of the 


voltage injected to the first phase of the 
second circuit. 
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Denoting the elements of matrices by double subscripted 
lower case letters and the matrix [Z]~^[i/j] by [U], forcing of 
the zero current conditions on the two unenergised phases of 
the second circuit, make equation (7.11) take the form 


u(5,l) u(5,2) ... u(5,6) 

u(6,l) u(6,2) ... u(6,6) 


M, 

M, 


M. 


0 

0 


(7.15) 


The open receiving end condition gives rise to another 
equation 

M = -[Tanh^A ] |L. (7.14) 

Replacing M in (7»13) by (7.14), we get 


u(5,l) u(5,2) . . . u(5,6) 

u(6,l) u(6,2) ... u(6,6) 


[TanhipS.] K 


(7.15) 


Denoting 


u(5,l) u(5,2) ... u(5,6) 

u(6,l) u(6,2) ... u(6,6) 


[Tanhil^Jl] by [UT] ‘ 


and substituting for known components of K by equation (7.12), 
(7.15) becomes 



0 



ut(l,l) ut(l,2) ..w ut(l,6) 




'd 

ut(2,l) ut(2,2) ... ut(2,6) 


^4 1 


_o 






( 7 . 16 ) 
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■where 


0 


0 

0 

0 


From (7.16), and K.^ can he expressed in terms of known 
quantities as follows: 




ut(l,5) 

ut (1,6 ) 

; -1 

ut(l,4) 

% 


nt(2,5) 

ut ( 2 , 6 ) 


ut(2,4> 


(7.17) 


and are jiothing hut the transforms of the induced 
voltages at the sending end on the second and third phases 
of the second circuit. 


Thesh transforms of the induced voltages on the 
floating conductors are to he transferred to time domain 
hy performing the n'umerical integration. For continuing 
the calculations, when the second and the third poles get 
closed, a n-umerically evaluated modified Fourier transform 
of these induced voltages on these floating conductors 
during the interval when they open shoiHd he added to 
the analytically known transforms of the second and third 
phase injected voltages of the second circuit, which come 
into the calculations from instants and respectively. 

At the instant the second pole gets closed and 

now the conditions to he forced at the sending end are known 
voltages on first and second phase conductors and zero 
current on the third unenergised phase conductor of the 
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second circuit, in addition to zero voltage conditions on 
the three phases of the first circuit. 


The known transformed voltages at the sending end 
determine the first five components of the vector K now as 
follows: 


= 0 ; K^ = 0 I = 0 j 

% = ®S2 + ^^32 


( 7 . 18 ) 


where Eg 2 1^ the analytically loaown modified Fourier 
transform of the injected voltage on the second phase of 
the second circuit (which starts from instant T^) and 7^2 
is the niimerically evaluated modified Fourier transform of 
the induced voltage on the second phase of the second 
circuit in the intei*val 0< t <T 2 ^. 

The zero current condition on the third phase of 
the second circuit gives rise to the ec^uation 


[ut(l,l) ut(l,2) ... ut(l,6)] 


0 

®sl 

(®s2-*-^s2) 


= 0 ( 7 . 19 ) 


where [UT] is the same as used earlier and 


0 = 


0 
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Prom equation (7.19), the unknown can be expressed 
in terms of known quantities as follows: 

Kg = -[ut(l,4) Eg 3 _ + ut(l,5)(Bg2 + ^^ 3 ) ]/ut(l,6) (7.20) 

Kg so obtained is the transform of the induced 
voltage on the third phase of the second circuit, valid in 
the time interval 0 <t addition of the numerically 

evaluated modified Fourier transform of this induced voltage 
in this interval to the analytically known transform of the 
third phase injected voltage is essential for carrying on 
the transient calculations after the instant T 2 » 

At instant T 2 » third pole closes and now the 

transformed voltages on all the phases at the sending end, 
valid in the entire time range are known. 

Hence for the open receiving end condition, an 
expression for the transformed voltage at the receiving 
end can be derived using the general expressions (7.10) 
and (7.11) as follows; 

= [SechiljA] 

where is the known transformed voltage vector at the 
sending end, 

A numerical integration of this transformed voltage 
to time domain gives the desired transient response. 
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This transient response, when added to the steady 
state response, calculated assuming the second circuit to 
he unenergised, yields the total response. 

7.7 COMPUTATIONAL RESULTS 

7.7.1 Double Circuit Versus Single Circuit Lines: 

The same double circuit line example, considered in 
Chapters 5 and 6 is taken up and the transient overvoltage 
at the open receiving end, resulting from the sequential 
closing of circuit breaker poles is evaluated by the procedure 
described in the last section. Next the presence of the first 
circuit is ignored and the receiving end transient voltage, 
resulting due to the injection of the steady state voltage, 
calculated at the sending end of the first circuit in the 
steady state analysis, is computed. The line is assumed to 
be transposed and the closing instants of the phases are 
taken to be 0, 3.3 and 6,7 m.secs. for both the cases. The 
phase sequence of breaker pole closing is assumed to be 
R’B'Y'. The 50 Hz parameters of the line are used in the 
calculations. The open receiving end voltage waveforms for 
the single and double circuit representations are compared 
in Figure 7.6, A maximum peak of 2,505 p.u. is observed to 
occur on phase Y’ at 8,3 m.secs. for the case of double 
circuit representation; whereas the maximm peak for the case 
of single circuit representation is -2.9 p.u, on phase B* at 

4.8 m.secs. Noticing the wide variation in the voltage 
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■waveforms and in "tlie maximum peak values (abouf 16 peicenf 

of bo'fch. bhe cases, if may be concluded fhaf flie 
double circuit modelling is essential for the energisation 
transient study. 

7.7.2 Transposed Versus Uhtransposed Double Circuit lines: 

Next the energisation transient of the same example 
for the untransposed configuration of the line is evaluated. 
The eigenvalue, eigenvector analysis at each frequency in 
the integration range is performed for the double circuit 
line to arrive at the modal transfozmiation matrix for using 
in these calculations ^ The voltage waveforms at the open 
receiving end for the transposed and untransposed configura- 
tions of the line are compared in Figure 7.7, Considerable 
difference in waveshapes is observed for the two configura- 
tions as there are six distinct modes and velocities of 
propagation for the case of 'untransposed line, whereas for 
the transposed line there are only three such distinct 
modes and velocities of propagation. A maxim'um peak 
voltage of 2,505 p.u, is observed to occur on phase Y' at 

8.3 m, secs, for the transposed line. The maxim'um peak for 
the untransposed line is -2,46 p,u, and it occurs on phase 
B* at 4.95 m.secs, 

7.7.3 Effect of Frequency Dependence of line Parameters: 

The effect of continuous variation of line parameter 

values with frequency on the transient overvoltage waveforms 
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due to sequential closing of transposed and untransposed 
double circuit lines is next examined. The line parameters, 
calculated at each frequency using Carson's formulae are 
used in the calculations. The resulting transient voltage 
waveforms at the open receiving end for the same example, 
with and without frequency dependence of line parameters 
incorporated for transposed and untransposed lines are 
compared in Figures 7.8 and 7.9 respectively. Again the 
closing instants of the phases are assumed to be 0, 3.3, 

6.7 m.secs. The phase sequence of closing is taken as 
R'B'y. For the transposed configuration of the line, 
without considering the frequency dependence of line 
parameters, the maximum peak overvoltage (2,505 p.u.) 
occurs on phase Y’ at 8,3 m.secs. When the frequency 
dependence of line parameters is considered the maximum 
peak value reduces to 2.21 p.u, (about 14 percent reduc- 
tion) and this occurs at 8,1 m.secs. on phase Y' . The 
corresponding peak values for the untransposed configura- 
tion of the line without and with frequency dependence of 
line parameters, incorporated, are -2,46 p.u. (on phase B* 
at 4.95 m.secs.) and -2.28 p.u, (on phase B' at 4.95 m.secs.) 
respectively. There is also a smoothening of oscillations 
in the waveforms when frequency dependence of line parameters 
is included for both the configurations of the line. 
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7.8 CONCLUSIONS 

In this chapter, a program to compute the energisa- 
tion transient arising on a transmission line, due to 
sequential closing of the circuit breaker poles, by the 
modified Fourier transform method, has been described. The 
method, first developed by Wedepohl and Mohamed^ using 
numerically evaluated Fourier transform of induced voltages 
on the unenergised phases in the calculations, has been 
extended here for the double circuit line case, A typical 
400 EV double circuit line example has been studied for 
energisation overvoltages for both transposed and untrans- 
posed configurations of the line. The need for the 
modelling of a double circuit line for energisation 
transient study has been established. The effect of 
frequency dependence of line parameters on energisation 
overvoltages, produced on a double circuit line, has been 
investigated. The introduction of frequency dependence 
of line parameters reduces the peak overvoltage and 
smoothens the oscillations in the voltage waveforms. 



CHAPTER 8 


SUMMARY ARD CORCLUSIORS 


8.1 C0RCLUSI0R5 

8.1,1 ; The study of electrical transients on EHV 
power systems has gained great importance » particularly in 
view of the increasing operating voltage of the transmission 
systems. The cost of insulation foims a major portion of 
the total cost of a transmission system. Especially as the 
system voltage increases to KJV level, the insulai:ion cost 
of the line increases at a very rapid rate. While doing 
these planning studies, if the models of the transmission 
line and other equipments that are en5>loyed in the transient 
calculations are not very accurate, it is necessary to use 
a high factor of safety to take care of the likelihood of 
the consequent margin of error in the calculated values. 

This ultimately reflects in higher cost of insulation. If 
this factor of safety and insulation costs are to be reduced, 
it can be done only if the calculated results predict the 
system performance with a high degree of accuracy. To achieve 
this, it is necessary to employ accurate models for all the 
components of a By#tem while doing the required studies. 

Bouble circuit lines are common features of any 
power system. Howevss™, of transients ■ 

on double circuit lines have not been reported. Nevertheless 
a few TNA studies fori, double circmt lines have been repeated. 
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as acknowledged in Section 5«1» An accurate mathematical 
model of a double circiiit line suitable for studying the 
electrical transients arising on it has been proposed in 
this thesis. Further, this model is made use of to stuiy 
the energisation and fault transients on a typical 400 KV 
double circuit line. This illustrates the use of the above 
method and also shows that it is essential to represent 
the double circuit lines adequately in any transient 
studies. 

There are various methods of calculating the over- 
voltage transients that are available in the literattnre. 
These methods and their special features are discussed in 
Chapter 1. Among these, the* modified Fourier transform 
method is regarded as one of the more acctirate methods, 
because it considers all the line parameters to be distri^ 
buted, can incorporate the effect of frequency dependence 
of line parameters with ease and can also handle 
untransposed lines. Therefore the proposed models in this 
thesis largely refers to the modified Fourier transform 
method. In one of the chapters, it has also been applied 
to Uram and Miller's method. Similarly, it can be extended 
to other methods, if desired. 

8.1.2: A mathematical model of a double circuit 
line, suitable for transiaat overvoltage studies has been 
developed in Chapter 5* It has been shown in Section 3*2.1 
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that there exists a constant modal transfoimation matrix 
for the case of a transposed double circuit line, vrhose 
conductors are situated at a mirror symmetric position 
with respect to a vertical plane. This transformation 
holds good at each frequency in the numerical integration 
range. Thus, for calculations of switching overvoltages 
involving transposed double circuit lines considerable 
saving in computation time can be achieved. On the other 
hand, in the case of untransposed double circuit lines, 
the eigenvalue, eigenvector analyses of two 3x3 matrices are 
to be performed, at each frequency, to arrive at the modal 
transformation matrix. 

8.1.3; The transient overvoltages, arising on the 
sound phases of a double circuit line due to the initiation 
of a single line to grotnad faialt on one of the circuits, 
have been studied in Chapter 4. Necessity for the study 
of the transients due to the initiation and clearing of 
fault in a system arises when effective methods are 
adopted for controlling the overvoltages, caused by 
energisation and re- energisation of the line. At present, 
as there are very effective methods, such as miiltistep 
resistance closing, controlling the closing instants of the 
individual poles of the circxiit breaker etc. available to 
limit the energisation and re-ehergisation overvoltages, the 
fault initiation and clearing transients may dominate the 
scene. These fault transients have been studied by various 
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authors in the past for single circuit lines using vario\xs 
digital methods. However » an accurate analysis of these 
fa\ilt transients on double circuit lines has not been 
reported. Moreover, in all the above studies, the lines 
were assumed to be transposed and the frequency dependence 
of line parameters have not been taken into account since 
they used lattice diagram or Laplace transform approach. 

In this thesis the modified Fourier transform method has 
been onployed for these studies. Further, this has also 
been extended to the double, circuit line model that has 
been proposed. A typical 400 EV double circuit line 
example has been considered. The transient voltage 
waveforms on the sound phases, following a sin^e line 
to groimd faiilt at mid span on the line, have been computed 
for both the transposed and untransposed configurations of 
the line. As shown in Figure 4.5, the maximum peak voltage 
in the case of untransposed line is observed to be nearly 
10 percent higher than that of the transposed case. 

In the above set up the presence of the second circuit 
is ignored for calculation of the same transients, thus 
representing a single circuit line. In this case fault is 
considered on one of the phases and the resulting transients 
on the remaining sp'^nd phases are evaluated. The voltage 
waveforms on the sound phases at the fault point for the 
double and single circuit representations have been shown 
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in Figure 4 •4* The voltage waveforms for the case of the 
double circuit line clearly exhibit more high frequency 
oscillations than those for the case of the single circuit 
line. This seems to be due to the fact that there are 
three distinct modes of propagation with three distinct 
velocities of propagation in the case of transposed double 
circuit line; whereas for ihe case of a transposed single 
circuit line, there are only two such distinct modes and 
velocities of propagation. A maximtim peak of -1.78 p.u. 
is observed to occur on phase-c’ of the second circuit 
for the double circuit line representation. For the 
single circuit line, the maximiuji peak of -1.98 p.u. 

(10 percent more) occurs on phase-c. This order of 
difference between the peak voltages for the single and 
double circuit representations shows that it is important 
to represent the double circuit line adequately. The 
effect of frequency dependence of line parameters on the 
fault transient overvoltages for both the transposed and 
untransposed configurations of the transmission line has 
also been investigated as shown in Section 4.6.3. The 
inclusion of th^ frequency dependence of line parameters 
is observed to reduce the maximum peak voltage by nearly 
11 percent and 13 percent (see Figures 4.6 and 4-7) for 
the cases of trajpasposed and •untransposed lines respectively. 
There is also a smoothening of oscillations in the transient 
voltage waveforms when frequency dependence is included in 
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the calcinations* The fault transients have been evaluated 
for various off-centre fault locations as shown in Table 4*2. 
The mid span fault is found to give rise to the maximum peak 
overvoltage. 

8 . 1.41 The modified Fourier transform method has been 
applied to analyse the energising transients on double circuit 
lines in Chapter 5* Algorithms have been developed to study 
the transients due to the closing of the first circuit with 
the second circuit unenergised, as well as for the transients 
due to the closing of the second circuit when the first 
circuit is operating in steady state. Here all the three 
poles of the circuit breaker are assumed to close simultane- 
ously during energisation. A typical 400 KV double circuit 
line, open at the receiving end has been considered and 
both the above mentioned energising transients at the 
receiving end have been calculated. The energising transient 
for the double circuit line representation has been compared 
with the transient caued by the energisation of the comes- | 
ponding single circuit line (ignoring the presence of the 
other circuit). For the untransposed configuration of the 
line, a difference of 20 percent in the peak transient 
voltages is observed for the single and double circuit 
representations (Figure 5.6). Apart from such a consider- 
able difference in peak overvoltages, the voltage waveforms 
for these two representations also differ widely in shape. 
This may be attributed to the reason that in the case of 
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untransposed double circuit lines there are six distinct 
modes and velocities of propagation; whereas there are only 
three such distinct modes and velocities of propagation in 
the case of untransposed single circuit lines. 

The effect of inclusion of the frequency dependence 
of line parameters is observed to be insignificant for the 
case of the simultaneous" closing of the breaker poles. This 
appears to be due to the absence of the zero sequence 
currents for simultaneous closing of the circuit breaker 
poles and only the zero sequence components of the resistance 
and inductance of the line (grcund return path taken into 
acco\mt) vary considerably with frequency. 

8.1.5 5 In. Chapter 6, the Uram and Miller's method 
has been extended to study the energising transients on 
double circuit lines. Uram and Miller's method is a time 
domain method, in which by making some assumptions, an 
incremental solution of the transmission line equations 
is made possible. The existence of a frequency independent 
modal transformation matrix is necessary for the application 
of this method. Such a transformation is shown to exist 
for the transposed double circtlit line in Section 3 ‘2.1. 

The same double circuit line example as in Chapter 5 has 
been considered. The transients at the open receiving 
end, following the energisation of the first circuit with 
the second circuit unenergised as well a^ that due to the 
closing of the second circuit when the first circuit is 
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operating in steady state have "been computed. The results 
obtained by this method are found to be closely matching with 
those obtained using the modified Fourier transform method. 
This serves as an important check on the validity and 
acc\iracy of the method proposed here as well a.s the two 
computer programs developed for this pvirpose. The Uram 
and Miller^s paethod takes much less computer time but can 
handle only transposed lines and cannot take into accotint 
the frequency dependence of line parameters directly. 

8.1.6: So far in all the studies all the three poles 
of a circuit breaker are assumed to operate simultaneously. 
However, in reality, there may be a difference of a few 
m.secs. in the operation of different poles of the circuit 
breaker. The incorporation of the nonsimixltaneous closing 
feature of the circuit breaker poles in the modified Fourier 
transform method of calculating energisation transients 
has been described in Chapter 7. The method has been 
extended to evaluate the transient arising due to sequential 
energisation of the second circuit of a double circuit line 
when the first circuit is operating in steady state. Th^ 
same 400 KV double circuit line example as of Chapters 5 
and 6 has been studied. When the sequential closing of the 
breaker poles is incorporated in the study, for the transposed 
configuration of the line, the peak voltage is 16 percent 
lower in the case of a double circuit line as compared to 
that of the corresponding single circuit line (Figure 7.6). 
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It may be noted here that there was no difference in the 
peak Values for transposed single and double circuit lines 
when the breaker poles are assumed to close simultaneously 
(see Section 5«6.2). The need for the modelling of the 
double circuit line has been established by comparing the 
transient arising on the second circuit obtained by ignoring 
the -presence of the first circuit with that obtained by 
the double circuit representation. 

The transient voltages at the receiving end 
following sequential closing for the transposed and un- 
transposed configurations of the double circuit lines 
have been studied. A considerable difference in wave- 
shapes has been observed for the two configurations as 
shown in Figure 7.7. 

The effect of inclusion of frequency dependence of 
line parameters on the transient overvoltages arising on 
the double ci^pcuit line following sequential closing has 
also been investigated. The introduction of the frequency 
dependence of line parameters is obsei*ved tO reduce peak 
overvoltages by 14 percent and 8 percent for the transposed 
and xmtranBposed configurations of the line respectively 
for the example considered. There is also a smoothening 
of oscillations in the voltage waveforms. 

,In this thesis, the transient arising due to re- 
energisation of a double circuit line on trapped charges 
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has not been considered. However, the method is general 
enough and this aspect may be studied in future. In the 
various examples considered in this thesis the source 
sides have been represented by lumped impedances. The 
values of these source impedances have been derived from 
the system short circuit levels. This simplified model 
for the source is somewhat approximate. However, here the 
aim is just to illustrate the application of the modified 
Pourier transform method for transient an^ysis on double 
circtiit lines and hence source sides have not been 
rigorously ropresented.- Howev'er, there is no 

limitation for the applicability of the method proposed 
here for a detailed source representaticm. The other 
lines terminating at the ends of the line tinder considera- 
tion, can be represented quite accurately by shunt elements 
at these buses, the values of which can be obtained by 
finding out the equivalent admittances offered by these 
lines at these buses, taking into consideration the distri- 
buted nature of these lines and the terminal conditions 
existing at the remote ends of these lines. 

In addition to the above the following suggestions 
are made for futirce extension of this work. 

1. Computation of transient arising due to the re- 
energisation of a double circuit line on trapped 
charge. 
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2. Fault clearing overvoltages and transient recovery 
voltages evalTiation for a system consisting of double 
circuit lines. 

3. Development of a general purpose program, using modified 
Fovirier transform method to accommodate general' soiirce 
side network. 

4. Studies of transients in systems, consisting of non- 
linearities such as transformer saturation, lightning 
arresters etc. using modified Fourier transform method 
■for the line. 

5*. Studies of overvoltage transients on double circuit 
lines with series capacitor compensation. 

S.2 SUMMAEY 

We have seen that in the case of a typical 400 E7 
line, considerable difference in calculated results is 
obtained when the line is represented by a single circuit 
or a double circuit model. In particular, the following 
results are very significant. 

1. The maximum peak overvoltage arising on sound phases 
following a single line to ground fault for the 
transposed double circuit representation is -1.78 p.u. 
The corresponding maximtim peak for the transposed 
single circuit representation is -1.98 p,u. 

2. For the transposed configuration of the double circuit 
line, maximum peak voltage due to f atilt initiation is 
-1.78 p.u. For the untransposed configuration, the 
corresponding maximum peak is -1.97 p.u. 
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3. Wlien freq^uency dependence of line parameters is consi 
dered, the maximum peak voltage due to fault initiation 
reduces by 11 percent and 13 percent for the transposed 
and untransposed configurations of the double circuit 
line respectively. ■ 

4. The maximum peak voltage due to simultaneous breaker 
pQ]Lg closing of the first circuit of a transposed 
double circuit line is -2.45 p.u. and occurs on phase-a 
at 9*2 m.secs. The corresponding maximum peak for 

the untransposed double circuit line is 2.24 p.u. and 
this occurs on phase— c at 12.8 m.secs. 

5. The maximum peak voltage due to simultaneous breaker 
pole closing of the second circuit of the transposed 
double circuit line is 2.02 p.u. and occurs on phase a* 

V at 1.02 m.secs. The corresponding maximum peak for 
the case of untransposed double circuit line is 
2*23 p.u. and occurs on phase a’ at 0.9 m.secs. 

6. With sequential closing of circuit oreaker poles 
incorporat ed > the maximum peak voltage due to second 
QpPQ-ait ©nergisation for transposed double circuit 
line is 2*505 p.u. and occurs on phase 1’ at 8.5 m.secs. 
The maximum peak value for the corresponding single 
circuit line is —2.9 p.u. and occurs on phase B’ at 

4.8 m.secs* 

7 . A maximum peak voltage of 2.505 p.u. is observed to 
occur on phase Y’ at 8.3 m.secs. for sequential second 
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circuit energisation of the transposed double circuit 
line} whereas the corresponding maxiniuirL peak for un— 
transposed double circuit line is -2.46 p.u. on phase B' 
at 4.95 m.secs. 

8. Reduction in peak values of 14 percent and 8 percent is 
observed in the transients due to second circuit 
sequential energisation for transposed and untransposed 
double circuit lines respectively, when frequency 
dependence of line parameters is considered. 

The order of magnitude of the error in the above 
cases is such that it appears vital to consider the followxng 

while doing any transient studies. 

a. Double circuit lines should be represented adequately. 

b. Transposition of the lines or otherwise, shotad be 
taken into account. 

c. Frequency dependence of line parameters can be 
ignored only when it is justified specifically. 

d. Non-simultaneous operation of the circuit breaker 
poles should be considered in these calculations. 

In the opinion of the author, optimum qosts of 
insulation of a transmission line can be achieved only when 
all the above factors are taken into account. These have 
been implemented in the computer programmes developed by 
the author and therefore the above suggestions are quite 
feasible. The validity of these suggestions is well borne 
out by the study of a typical 400 KV line that has been 



carried out throughout. In particular power syetens, the 
errors due to the various epproximations --Eiay he even higher 
than those reported here. Therefore, the importance of 
accurate representation of the components of the system in 
transient calculations cannot he imder emphasised. ^Methods 
to incorporate all the above suggestions have heen 
successfully implemented in the procedures proposed in this 


thesis. 
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APPENDIX A 
SYSTEM DATA 


length of the line 
Source inductance 
Earth conductivity 
50 Hz line parameter 


240 Kms. 

0.1 H/phase 
100 Ohm-m. 

matrices for the untransposed 


c onf i- 


gur at i on are as f oil ows : 



0.079658 

0.049180 

0.048136 

R = 

0.049180 

0.078826 

0.048053 


0.048136 

0.048053 

0.077874 


0,002493 

0.001132 

0.000829 

L = 

0.001132 

0. 002434 

0.000979 


0.000829 

0.000979 

0.002291 


1.89230 

-0.35843 

-0.10716 

C = 

-0.35843 

1.81560 

-0.34268 


-0.10716 

-0.34268 

1.74110 


X 0.62 ohm/Em. 


X 0.62 H/Km. 


X 0.62 X 10“® f/Kin. 


For the case of the tranposed line, the line parameters 
are obtained by finding the average values of the diagonal 
and off-diagonal elements of these natrices separately, 
corresponding to the untransposed case. 
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APPENDIX B 

IINE PARAMETERS CAlCUlATION - CARSON'S FORMULAE 


Self impedance of conductor = + 0.00159f 

+a0.004657f log 

&MR 

Mutual impedance between two conductors = 0.00159f 

+ j 0.004 65 7f lo^ 2M91L£JI1 

where = resistance of the conductor per mile 

f = frequency in Hz, 

P = earth resistivity in ohm-m 

GMR = geometric mean radius .of the conductor 
in feet 

<iab “ distance between the conductors in feet 
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APPENDIX G 

SYSTEM data - BOONYUBOL‘S EZAMPIE 

Positive sequence resistance of the line = O.O 3 ohm/mile 
Positive sequence inductance of the lin6= 1.56 mH/mile 
Positive sequence capacitance of the line = 0.0189 lif/mile 
Zero sequence resistance of the line = O. 5 I 6 ohm/mile 

Zero sequence inductance of the line =5.76 mH/mile 

% 

Zero sequence capacitance of the line = 0.0138 |Af/mile 
Frequency of supply = 60 Hz 

length of the line = 180 miles 

The fault is simiilated at mid span. 

The source impedance has "been taken as zero. 

The fault has been assumed to occur on phase 'a' at the 
instant wHen phase ’a‘ voltage is at positive peak. 
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APPEm)IZ D 

SYSTEM DATA PCR PAUIT STUDIES 


LengtE of the line 480 Ems. 

The line is fed on both the ends hy similar sources, which 
are purely inductive of value 0.05 H. 

Earth conductivity = 100 ohm-m. 


The computed line parameter matrices at 50 Hz for 


the untransposed configuration are given below. As these 
matrices are symmetric, only the lower triangular elements 


are given. 


[R]= 


0.079658 

0.049180 0.078826 
0.048136 0.048053 
0.049658 0.049180 
0.049180 0.048826 
0.048136 0.048053 


0.077874 

0.043136 

0,048053 

0.047874 


0.079658 

0.049180 0.078826 
0.048136 0.048053 0.077874 


X 0.62 ohm/Km. 



0.002493 

0.001132 

0.000829 

0.000972 

0.000876 

0.000754 


0.002434 

0.000979 0.002291 
0.000876 0.000754 
0.000848 0.000786 
0.000786 0.000829 


0.002493 

0.001132 0.002434 
0.000829 0.000979 0.002291 


X 0.62 H/Km 
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1 . 39230 

- 0.35843 1.81560 

- 0.10716 - 0.34268 1.74110 

- 0.12082 - 0.08571 - 0.06246 1.89230 

- 0.08571 - 0.10012 - 0.12187 - 0.35843 1.81560 

- 0.06246 - 0.12187 - 0.26750 - 0.10716 - 0.34268 1.74110 


X 0. 62 X 10 ® f/Km. 


For the transposed configuration, the self elements 
of the line parameter matrices are obtained by taking the 
average value of the self elements of the corresponding 
parameter matrices of the untransposed case, given abo-ve. 
Similarly, the mutual elements between the phases of 
each circuit and the mutual elements between the two 
circuits for the transposed line are the average values 
of the corresponding elements of the untransposed case. 
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APPENDIX E 

UHAM AKD MIILER'S METHOD - THEORY 


The Uram and Miller's method is a time domain method^ 
in which by making certain assumptions, an incremental 
solution of the transmission line eq^uakcions for the voltage 
and current at each instant at any point on a transmission 
line, subsequent to a switching operation is made possible. 


The voltage and current vectors at any point on a 
transmission line, subsequent to a switching operation are 
governed by the equations 


3 2(x»t ) 
” l“x 

3 i(x,t ) 

sTx 

neglecting the 
line. 


3 i (x,t ) 

= [R] i (x,t) + [l] -j-i 


= [C] 


3Z(^»‘t) 


3t 


conductance to ground of the transmission 


Taking the laplace transform with respect to t» 
these equations become 
dV(z, s) 


dx 


= { [R] + s[Lj}I(x,s) 


(E.l) 


dl(x,s) 

dx 


s[C] Y(x,s) 


(E.2) 


These two equations can be combined to get 


d^V(x,s) r Tr T / N 
— [z]Cy] I(x,s) 

dx^ 


(E.3) 
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whore [z] ={ [r] + s[Ij] } 
and [y] = s[cj 

The above coupled matrix differential eqixation 
can be decoupled by performing the modal transformation, 

V(x,s) = [t] F(XfS) 

where [T] is an eigenvector matrix of the matrix product 
[Z][Y]. 

Equation (E.3)» after the modal transformation, takes 
the form 

d^F(x,s) n 

= [T]~^ [z][y][t] E(x,s) 

dx^ 

where [t]‘"^[z][y][t] is a diagonal matrix. 


The solution of this equation for the ith modal 
component can be writt en as 


F^(x,s) = 


Kli(s) e 


-f(Z, 


mi 




K2i(s) 




Y . )x 
mi 


,i-l,lT. 


where and are constants, independent of x 

and Z . Y . = ith diagonal element of the diagonal 
mi mi 

matrix [T]-^[z][T][t] . 

In the vector notation, the solution can be writt ©a as 
F(x,s) = K^(s) + Kg,(s) 


where 
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e; = [K^^(s) e ’ ... ^ 


and 


V(Z„n 7„t)^ 


„ , X + rtr r ^ ""ml " 111'"' „ , 1 

1^2 (s) e^ = LK:23 _(s) e ... £211^®^ ® J 


Transforming back to the phase quantities 
V(z,s) = [t] 1(z,s) = [T] {K. (s)e“ + K (s)ej } 


Prom equation (E.l), the solution for the transformed 
phase current vector can be written as 


I(x,s) = [t] [ a]"^ { K^(s)e“ - E (s)ej} 


+7 


where [^] = 


0 ... 0 

0 «Zj,2/T^2) ... 0 

«• ««•«« • ••••• 

0 0 ... f(z^A,a) 


Prom these general transformed voltage and current 
vector solutions, the expressions for the sending end and 
receiving end transformed voltage and current vectors 


can be written as 

Ig(s) = [t] {E^(s) -fK^Cs)} (E. 4 ) 

Ig(s) = [T][fi]“^ {K^(s) - 12 ( 3 )} (E. 5 ) 

I^{s) = [t] { 4 (s) + B(s)} (E. 6 ) 

XL 

Ig(s) = [T][n]“^{ l(s) - fiCs)} (E.7) 
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where the ith components of the vectors A(s) and £2^®^ 
related to the ith components of K^(s) and .B(s) by the 
expressions 


( s ) — ( s ) e 


-f(Z . Y .) Z 
* ' mi mi 


-f(Z . Y .) li 

J rr / ^ -O ^ ^ mi mi 

and E^i ^ ® B ( s ; e 

In the above expressions 


a = the length of the line 

= e«\i C^)(l 

mi 

where R_ - » L . and C . are the resistance, inductance and 
mi mi mi 

capacitance per ■unit length respectively offered by the 
line for the ith mode of propagation. 

If ““I ’ "■ ^ '^‘4’ 


r(^) = f(^) (1 + 


R . -4- 

mi ' 2 


Z . 
mi 


mi 


iir. 

mi 


) 


If 

mi mi ml 




mi' mi 


Hence A>, (s) = K-, • (s ) e 

^ (B.8) 

K2i<s)= B^(s) e “f e 
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When equations (E.4-) - (ll»7) are transformed back 
to time domain, the sendiag end and receiving end voltage 
and current rectors become 


IgCi:) = [f ] ' (t ) + k^lt) } 

(E.IO) 

igCO = ^ 

(E, 11) 

LijjCt) = [T] {a(t) 4-;b(t)3 

(E.12) 

ijj(t9= {a(t) - bit)} 

(B.15) 


where [fi] = 

° 0 


0 O .... 1 


By the shifting property of Laplace transforms, 
when f'onctions and are related by expressions as 

in equations (E.8) and (E.9) to the functions 
Bj^(s) in the Laplace domain, the corresponding inverse 
transformed functions in time domain are related by the 
expressions 


/2)f(C ./I 

a^(t) = e ^ u(t-T^) , i=l,N 


k (t) = b O-t ) uC-t ) , 1=1,N 

2i j. -L X 
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where = the travel time of the ith mode of prippagation 
of the travelling wave to traverse the length 
of the line 

Thus the voltage and current vectors at each instant 
of time at any point on a transmission line, subsequent to 
a switching operation can be expressed in terms of the 
values of the four function vectors k^(t), k^Ct), a(t) and 
b(t) at that instant. Among these four function vectors, 
a(t) and k^Ct ) have been shown to be the attenuated delay 
functions of the other two function vectors k^(t) and 
b(t) respectively. The function vectors k^(t) and b(t) 
can be expressed in terms of the function vectors a(t) 
and k^Ct) by the terminal conditions of the line. For 
the case of the line, which was relaxed before the 
switching operation, initially till thed;ime equal to 
the travelling time of the line, the delayed functions 
are zero. 


i. e. 

II 

O 

i = 1,N 

for 

V 

o 

< T. 
1 

and 

k2i(t) = 0, 

i = 1,N 

for 

0 < t 

<Ti 


Thus an incremental solution of the transmission 
line equations in time doitaiii is possible for an initially 


relaxed line. 
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APPENDIX E 

PLOW CHART POR COMPUTATION OP SEQUENTIAL CLOSING TRANSIENTS 
BY MODIPIED POURIBR TRANSPORM METHOD 


Read Data 


a» At, Aw* 

[q]» [Q]”^* Line data 



1 

0) = (t> + Ao) 


YES 







nil 

= 

1 



0) 

= 

0 



T 

= 

T-, 



max 


1 



l<l 

01 

II 

0, 

» 0 

<t 


2g(t) = 

0. 

» 0 

< t 

< T_ 


&J = 0 
T„„„ = T, 


[Y (01)] = 

i/[Rg+(a+3 a))LgJ 0 0 
0 0 0 

0 

O' 

0 

J 


<t < T, 


/[Rg+(a+Ot*>)Lg] 0 


[Yg(w)] = 


l/[Rg+(a+ 3 w)I-g] 0 


[Yg(w)] = 


/[R„+(a+3w)I'„J 0 

S o 

0 l/[R+(a+3 01 ) 1 - 3 ) 
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APPENDIX G 

0BRA”1UCKN0W LINE DATA 


Length of the line 400 Kms 

Shunt reactors of reactance value 3200 ohms at hoth the ends* 
Source reactance at Ohra end is 67*8 ohms, which takes in- 
to acco-unt "both the local generation, short circuit power 
supglied through transmission lines on l.ower voltage levels 
and the short circuit power of the 400 EV transformers. 

The line parameters are as follows: 


Positive sequence resistance 0. 

Positive sequence inductive reactance 
Positive sequence capacitance 
Zero sequence resistance 
Zero sequence inductive reactance 
Zero sequence capacitance 


027 ohm/Zm 
0.331 Ohm/Km 
10*95x10“^ f/Em 
0. 261 ohm/Em 
1.31 ohm/Ein 
8.42x10”^ f/Em. 



181 


CURRICULUM VITAE 


1. Candida't e * s names R» Balasubramanian 

2. Academic Background; 


Degree Specialisation 

B.E. Electrical 

Engineering 

M. Sc. (Engg. ) Power 

Syst ems 

3* Publications; 

i) Calculation of transient due to fault initiation 
on a double circuit transmisaion line, Accepted 
for publication in the Proceedings lEB, London. 

ii) Calculation, of energising transient on a double 
circuit transmission line, Submitted for publication 

in Proceedings IBB— IEEE (India). 

iii) Accurate evaluation of energising transients on 
double circuit EHV tranwbission lines. Submitted 
for presentation in IEEE Winter Power Meeting, 1977. 


Institute 

Year 

P.S.G. college of 

1968 

Technology , Madras 


University 


Annamalai 

1970 

University 




i 



511S1 

This book Is to b@ returned on tho 
date last stamped. 















1 







r 

CD 4 72.9 

EE" G - D' 




